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Abstract 

In F-theory GUTs interactions between fields are typically localised at points of enhanced sym- 
metry in the internal dimensions implying that the coefficient of the associated operator can be 
studied using a local wavefunctions overlap calculation. Some F-theory SU{5) GUT theories 
may exhibit a maximum symmetry enhancement at a point to Eg, and in this case all the oper- 
ators of the theory can be associated to the same point. We take initial steps towards the study 
of operators in such theories. We calculate wavefunctions and their overlaps around a general 
point of enhancement and establish constraints on the local form of the fluxes. We then apply 
the general results to a simple model at a point of Eg enhancement and calculate some example 
operators such as Yukawa couplings and dimension-five couplings that can lead to proton decay. 
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1 Introduction 

Given that string theory, if correct, should account for all current and future particle physics 
experimental data, string phenomenology is a data-rich subject. However much of this data 
lies not with the observed states of the Standard Model and its extensions, but with the op- 
erators of the theory. This includes measured values of operators, or couplings, in the SM as 
well as constraints on operators of possible extensions of it to, for example, supersymmetric 
or Grand Unified Theories. Quantitative calculations of operator coefficients within realistic 
string theory models are typically difficult because they require tools beyond algebraic geom- 
etry which are difficult to apply to complicated higher dimensional geometries with various 
fluxes and branes present. This means that, apart from for particularly simple geometries, 
it is difficult to say much beyond the bottom-up rule that if an operator is not forbidden 



by some symmetry then it is present with an order one coefficient. This approach can be 
refined by considering for example the scaling of operators with various moduli and other 
coarse properties of the background, and for many applications this is sufficient, nonetheless 
a quantitative computation is often missing. 

One class of constructions where we might be able to present a quantitative study of 
operators within a general setting is in type IIB or F-theory where the set of cubic normalis- 
able operators can be calculated using wavefunction overlaps on smooth backgrounds [n-fT5]rl 
This by itself is still a task beyond our reach for a general background, but there is an ad- 
ditional advantage which is that the couplings are naturally associated to only a relatively 
small patch within the full Calabi-Yau. This follows for the common constructions where 
gauge fields are associated to divisors wrapped by 7-branes, matter fields to curves where 2 
branes intersect, and cubic couplings to points where 3 curves intersect |19fi25j . Given this 
we can study the couplings by calculating the wavefunctions and their overlaps on a local fiat 
patch around a point of interaction. This was done for Yukawa couplings in [ll[71[9l[TT] and 
for some other couplings, particularly to heavy modes, in |13j. 

There is no doubt that one of the most attractive aspects of IIB and F-theory models 
is the control over the full global setup. They are the constructions where the vacuum 
structure and its stability can be argued to be best understood [26], and there has been 
substantial development in understanding of global model building in F-theoryo However it 
is also fair to say that one of the deepest problems facing string phenomenology is that of 
the landscape and in particular how to extract semi-universal phenomenological aspects. A 
major advantage of the local approach in this respect is that the large parameter space of 
the CY geometry and the many fiuxes gets pulled back, or projected, to just a small set of 
adjustable local parameters around the point of interaction. This is particularly appealing in 
the sense that although the large number of vacuum configurations all necessarily affect the 
vacuum energy, or cosmological constant, the projection to the local patch implies that the 
various configurations can give rise to a comparably small set of particle physics models. We 
can then hope to use this and the additional calculability of local models to impose a large 
number of constraints, corresponding to the many operator coefficients, on a relatively small 
number of local parameters. 

An important disadvantage of studying the theory around a point of interaction is that 
the local fiuxes and geometry at different points, say at an up-type Yukawa point and a down- 
type Yukawa point, are not related in a simple way and so it is difficult to calculate relations 
between the two, such as the CKM matrix. It is therefore interesting to consider models 
where both the Yukawa points are located in the same local patch. Indeed if the observed 
fiavour structure of the SM, which has strong correlations between the up and down sector, 
has a geometric origin it is natural to expect that the two points should be experiencing the 
same geometry. Since the interaction points are associated to an enhanced local symmetry, 
for example 50(12) and Eq for S'C/(5)-GUT Yukawas, their coincidence corresponds to a 
further enhancement in the symmetry. For the Yukawas this amounts to at least Ej but if we 
consider more interactions it is natural to consider the maximum enhancement to Eg . Indeed 
this reasoning was the motivation presented for studying the idea of a point of Eq enhanced 
symmetry in [34j. Such a construction enjoys this motivation, but a further slightly different 



^The singular limit of geometric IIB/F-theory constructions which are branes on singularities also allow 
for calculating couplings but using different techniques |14fll8) . 

^See ^5] for reviews and [27H33) for some of the most recent developments. 



motivation is that in the spirit of being able to calculate operators on a local patch we are 
motivated to study the richest and most complete theories to which these calculations apply. 
This way we are able to calculate many operators that in models where the interactions are 
spread out throughout the CY we would not be able to. In some ways these are the simplest 
realistic theories where we can hope to calculate the bulk of the operators as a function of a 
relatively small set of parameters. 

The work in this paper is motivated by the idea of a highly calculable, quite general, 
theory around a point of Eg- The idea being to take the point of £"8, or some other enhanced 
symmetry, as a starting point and calculate the couplings of the fields to each other as a 
function of the local values of fluxes and geometry. There are GUT singlets typically present 
at such points which parameterise deformations of it, and therefore these theories also include 
the set of theories where the singlets have a relatively small vev, or where the point of E^ is 
slightly deformed. In principle such theories are rich enough to have substantial applications 
to phenomenologyo 

Aside from applications to models based on the point of Eg we will study quite general 
properties of the effective gauge theory and its wavefunction solutions. Indeed the bulk of 
the paper will not specify a particular gauge group but solve for the most general local forms 
of the wavefunctions for arbitrary gauge flux and Higgs backgrounds, in the leading linear 
approximation. Once these are obtained we will apply them to a toy model around a point 
of -Eg. The paper is set out as follows. In section [2] we will set up the general formalism that 
we will use. We begin by introducing the effective theory and discuss in detail constraints 
required for its validity. We then study quite generally the relation between the local flux 
and Higgs values, and the form of the wavefunctions and what type of chiral matter they 
correspond to, thereby developing the local relation between flux and chirality which plays 
an important role subsequently. In section [3] we present explicit solutions for the form of 
the wavefunctions, for both massless and massive modes, and their overlaps for quite general 
flux and Higgs backgrounds. In section [4] we utilise the results developed in the previous 
sections to study a toy model based on a point of Eg enhancement. Within this model 
we calculate some phenomenologically interesting quantities like the Yukawa couplings, and 
various higher dimension operators particularly those associated to dimension-five proton 
decay. We summarise the work in section [5l In the appendix we present a semi-local analysis 
of the flux and chiral spectrum in the absence of any monodromies, so for a fully split spectral 
cover. 



2 The effective theory 

The setup we are studying is of a 7-brane in F-theory, wrapping a four-cycle S, carrying a 
gauge group Gs and exhibiting an enhancement at a point on its worldvolume to a gauge 
group Gp which is at least 2 ranks higher. In the upcoming sections we will keep all the 
formalism general so that it can be applied to any 7-brane gauge group enhancing up to any 
higher group. However for concreteness we will refer to the particular case of Gs = SU{5) 



■^However there are two effects which may, but not must, play an important role in phenomenology which 



we do not consider, these are local monodromies [51I35J and possibly non-commutative deformation associated 
with flux or non-perturbative effects ^,'9l,'36lirT|. Calculating the operators in a setup including these effects 
most likely will require numerical analysis which we leave for future work and in this paper we will consider 
as a starting point the most basic scenario where these effects are absent. 



enhancing to Gp = £"§, and in section[3]we will consider models that are specific to this choice. 
In the infrared the theory on the 7-brane is described by a twisted 8-dimensional M = 1 gauge 
theory, with gauge group G and support on M^'^ x S, where S is a 4-diniensional Kahler sub- 
manifold of the F-theory Calabi-Yau fourfold X. The twisting accounts for the embedding 
of the 7-brane into the four-fold in a supersymmetric way |2HI22j. The field content of the 
theory can be written in terms of 4-dimensional A^ = 1 multiplets as 

*^mn ^ \\''fH)mn ^Xmm limn] ; 

V = {7j,A^,V} , 

where by abuse of notation we denote here and henceforth both the 8-dimensional fields and 
their internal wavefunction profiles by the same symbols and will only make an explicit dis- 
tinction when necessary. The subindices on the fields denote their local differential structure 
on S and so take values m = 1,2. Here A and ^ are chiral multiplets with respective F-terms 
G and Ti. V is a vector multiplet with D-term V. A^ and i'^H)mn ^^^ complex scalars while 
i/jrh, Xmn, and 7] are fermions. 

Following the ideas presented in |20p22) we can model the intersecting brane combinations 
that lead to an enhancement to Eg at a point in S by considering the 8-dimensional gauge 
theory to have the gauge group G = Eg. Within this theory we give the scalar Higgs 
component {<fH)mn ^ spatially varying vacuum expectation value (vev) which, at a generic 
point on S breaks the gauge group to SU{5)gut x U{1)^. The Higgs vev vanishes on certain 
curves and points in S on which therefore the gauge group enhances. Relevant to us will 
be curves where the Higgs vev vanishes along one U{1) factor inducing an enhancement of 
SU{5) by rank 1 to SU{6) or SO (10). As well as a point where the full Higgs vev vanishes 
inducing an enhancement by rank 4 to Eg. 

The equations of motion of the theory on M.^''^ x S were derived in [211122^ 19]. Solving 
them on some complicated background geometry S is a prohibitively difficult task. However 
since our primary interest is in the region around a single point, of Eg enhancement, in S 
we can approximate the theory locally as flat space S" ~ C2. Within this approximation the 
equations simplify considerably. 

Having described the general setup we go on to specify the details of the model. These are 
similar to, but more general than, the setup discussed in |T3| and the related configurations 
of [37H39t [TT]. to which we refer for further details on some of the concepts introduced below. 

2.1 The equations of motion 

We are interested in the equations of motion for fiuctuations of the fields about a background 
with a Higgs vev and non-vanishing gauge field flux along S. Locally, the leading order 
contributions from the Higgs and gauge fields are linear in the local co-ordinates zi and Z2 
and so take the form 

{ipu) = MKRm^Zi Qadzi A dz2 + . . . , (2.2) 

{A) = -Mk lui{Mt,z,dzj)Qa + . . . , (2.3) 

where the dots denote higher order terms in the two local complex coordinates zi and Z2- 
Here the mf and M?- are complex constants which, in the case of £^8 enhancement, denote the 



vevs along the 4 L'^(l) factors with associated generators Qa- Actually it is more convenient 
to turn on the vev within S{U{1)^) so that a = 1, ..,5 but we have to impose an additional 
tracelessness constraint 

5 5 

Y^n^ = Y,Mt^=Q. (2.4) 

a=l a=l 

We will also allow for flux MJ- along the Hypercharge U{1) inside SU{^)gut with generator 
Qy- Note that the local expansion of the Higgs and the flux begins with a linear term in 
the Zi and there is no constant term. For the Higgs background this amounts to defining the 
enhancement point to be at the origin. For the gauge field a constant term locally is pure 
gauge and so can be gauged away. 

The mass scale Mk is a local mass scale which involves the cutoff scale of the theory M^, 
and scales with a homogeneous rescaling of the local metric by a length scale R\\ as 

Mk = ^. (2.5) 

II 

The dimensionless scale R is associated to the local normal length scale R± to S and scales 
as 

R = R\\Ri_ . (2.6) 

The scaling of the Higgs with R± follows from the pullback of the normal metric component 
to the world- volume of the 7-brane as discussed in [13] . The gauge field is real and leads to 
flux of complex type (1, 1) and the Higgs is holomorphic, these two requirements amounting 
to the preservation of half the supersymmetry by each one. Further, given the local flat 
metric with associated Kahler form 

to = - {dzi A dzi + dz2 A dzg) j (2-7) 

the D-term constraint, uj A F = 0, imposes that 

Re Mfi = -Re M^^ ■ (2-8) 

Note that the flux associated to the gauge field (|2.3p reads 

F = -M|. {Mij + liiji) dzi A dzj . (2.9) 

Not all the flux parameters appearing in (j2.3p also appear in (j2.9p : for example only the real 
part of Mil and M22 appear in the flux. However the wavefunction calculation does depend 
on the choice of gauge field and so these parameters do feature in some of the calculation. 

With this background the equations of motion for fermionic fluctuations which lead to 
massless four-dimensional fields read [2T | [22 | [9l [39lll 11113) i 

D"^ = 0, (2.10) 



Note that we have sHghtly changed conventions from J^ to match those of [TT] by flipping the sign of 
some components in D . 
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(2.11) 



and 



D7 ^d,- -qaM^^-z, 
D7 = -RqafhiZi 



D+ = a, + -QaM-^z, 



1,2 



^3+ 



Rqamfzi 



(2.12) 
(2.13) 



Here qa denote the charges of the fields under the generator Qa- Since the charges are always 
contracted with the Higgs or fluxes it is convenient to introduce the notation 



M, 



^J 



qaM^, , rui = qaTJii . 



«J 



(2.14) 



It is useful to consider the 8-dimensional fields as a vector ^ which contains the ferniionic 
components of a single A/" = 4 niultiplet. The reason is that four-dimensional mass eigenstates 
will in general be combinations of the 8-dimensional fields. The particular splitting (j2.ip 
is rather more convenient for writing down the equations in coordinates. The equations 
of motion for fluctuations which lead to massive four-dimensional fields can be written as 



"\I/ = Inixf^ 



(2.15) 



with associated four-dimensional mass M^mx- 

Having described the effective theory we discuss the range of parameters for which the 
theory is valid. This was studied in |13] which we summarise and expand upon. There are 
a number of approximations being used and it is worth splitting the associated constraints. 
Firstly we have the constraint that the metric takes the flat form (|2.7p . Note that with our 
conventions we first write it as a flat space metric and then rescale away the factors of ii!|| by 
writing in terms of a vielbein basis. This means that corrections to the metric are suppressed 
within the region defining our local patch 



Patch : \zi\ < 1 



(2.16) 



Outside this patch we can not trust the form of the equations and their associated wavefunc- 
tion solutions. 

Even within the local patch we may have corrections to the theory itself. These come 
from large localised energy densities which are associated to derivatives of the Higgs and 
background gauge field and so to keep these small we require 



M 



«j 



m 



<i 



miRj_ 
Rii 



<. 1 



(2.17) 



These amount to dilute fluxes and small brane intersection angles and if they are violated we 
expect corrections to the theory coming from higher derivative operators. 

The full constraints are most simply satisfied by taking 



i?ll » i2_L > 1 . 



(2.18) 



As discussed, it is important to keep in mind that all of the parameters, -Rii, R±, rrii and 
Mij are local parameters and so are not directly tied to their associated global scales. Even 
so, at least in the local patch the manifold should be near homogeneous, and if this near 
homogeneity is extended to the whole of Squt and the CY then we have that 



1 



1 



R\\ ~ \o\{SGUTy^ - a^tjT ^ (2-19) 

R^ ^ {^^^Y . (2.20) 

where p is some power which varies according to the embedding of Squt hito the CY. For 
example p = 1 for a torus and 1/3 for a contractible Squt- 

Given this there are a number of important observations. The first is that the decoupling 
limit: Mpianck -^ oo with acuT -^ finite, is not naturally compatible with the constraint 
coming from small Higgs variations (|2.17p . This is because the parameters mj, like the fluxes 
Mij, are quantised in a global context. Physically this issue amounts to the fact that if the 
C/(l) branes go into the bulk of the CY they can not have small intersection angles with the 
GUT brane throughout the surface Squt- Intuitively this can be seen by noting that the 
Higgs field measures the distance to the bulk branes and so in going around a loop in Sgut, 
of typical length Ru, it must range all the way to the global scale R± and back again and 
therefore must have a large gradient and so large intersection angles. In this case the approach 
of modeling the intersecting brane setup using an 8-dimensional theory with a higher gauge 
group which is broken by the Higgs field breaks down. Rather the appropriate model would 
have to be 10-dimensional with a product gauge group treating the bulk branes as separateO 

However it is important to stress that this does not imply a breakdown of the theory for 
the actual finite value of Mpianck- Rather it requires some degree of inhomogeneity in the 
CY so that the local R± in patch differs enough from the global value which gravity feels. 
In fact this is another motivation for considering the point of Eg since then all the relevant 
calculations depend on the local value of R±. Another possibility is that the intersecting 
brane configuration is such that the U{1) bulk branes also wrap local contractible cycles so 
that the Higgs does not range over the full CY volume. 

Another consequence of extending the local scale R\\ to Sgut is that the local approx- 
imation, where the wavefunctions are localised on a small patch within Sgut restricts the 
parameter space significantly. This is because of the relatively strong expected coupling at 

the GUT scale so that Ru ~ ctGUT ~ ^ (^) ^^ J^^^ ^°^ that large (note that additional mas- 
sive fields which typically increase the GUT coupling make this situation even worse). To 
localise wavefunctions within patches much smaller than this requires strong energy densi- 
ties for which our effective theory would break down. This means that the wavefunctions 
sometimes still have a significant tail at distances \zi\ ~ 1 from the point of Eg and so for 
those wavefunctions the accuracy of the overlap calculation is reduced. Practically however 
these corrections can be kept well under control. On the other hand, as discussed in the 
introduction, the smallness of Sgut implies that even for general F-theory models which do 
not have a point of £"8 enhancement, since the wavefunctions on Sgut may have significant 
overlap, a model based on a point of Eg may be a decent approximation. In the rest of the 



''A more field-theory perspective on tiiis issue can be seen by tliinking about tiie energy scale at which the 
theory enhances to the higher gauge group all over Sgut- This energy scale goes like the bulk scale, for a 
non-compact setup diverges, and in particular is higher than the local cutoff scale of the theory and so such 
a theory based on an enhanced group is inconsistent. 



analysis we can consider working in the acuT — ^ limit where the results discussed and the 
localisation (onto matter curves) approximations become exact. 

Just as the Higgs field induces localisation onto the matter curves the flux on the matter 
curves induces localisation within the curves. However this latter effect does not become 
exact in the limit acuT -^ but rather the degree of localisation is controlled directly by 
the flux M. Since this is expected to be of order one this localisation is typically less strong 
than that transverse to the matter curve and forms again a source of possible corrections 
to the theory coming from the wavefunctions sensing the curvature corrections of the metric 
on the matter curve. The localisation along the curves is not so important for calculating 
the operators of the theory since the appropriate integrals are fully localised within Scut 
by the Higgs vev. However it does feature in the normalisation of the wavefunctions and in 
determining the local chiral spectrum of the theory. We discuss this in more detail in the 
next section. 

2.2 Local and global modes 

The key objects in calculating the operators of the four-dimensional theory are the internal 
wavefunctions of the four-dimensional fields along Sgut- The background Higgs profile ()2.2p 
ensures that the wavefunctions are localised onto matter curves. The profiles of the wave- 
functions along the curves themselves are more complicated and are primarily sensitive to 
the background flux (j2.3p . As well as controlling the wavefunction profiles along the matter 
curves the fiux also generates chirality in the massless spectrum. The interplay between these 
two effects of the flux plays an important role in the spectrum analysis and particularly so 
for the case of a point of Eg . In this section we study and develop this connection as well as 
elucidating more generally the question of how much can we learn about the wavefunctions 
from a given local patch along the matter curve. 

2.2.1 A global example 

As a reference to the general concepts we will discuss it is useful to have an explicit example. 
In this subsection we briefly study some properties of wavefunctions over an explicit matter 
curve C which we take to be complex projective space P^. This model is slightly different 
from our model in that it is 1-complex dimensional and untwisted, where as in our model 
wavefunctions are localised onto 1-complex dimensional curves inside a 2-complex dimensional 
twisted theory. Nonetheless the essential features will be similar. 

We begin by briefly recalling the most common way to count chirality on P^ with back- 
ground flux. We specify the background flux through the transition functions of the line 
bundle which for homogeneous coordinates Zi, and integer flux M, go as (^Zf/Z^^ , with 

patches labeled by a and b. Therefore global holomorphic sections are counted by homo- 
geneous polynomials of degree M — 1 of which there are M. Then, since we are counting 
solutions to the F-term equations only, we can go to the holomorphic gauge for which the 
Dirac operator reduces to the simple Dolbeault operator and so solutions are given by any 
holomorphic sections. There are M of these and so there are M massless modes, which are 
chiral since the conjugate partners would correspond to anti-holomorphic sections of which 
there are none because of the holomorphic line bundle transition functions. 



This way of counting chirality is intrinsic to the holomorphic gauge. However the holo- 
morphic gauge is not physical and is only useful for understanding rather coarse properties. 
For example it does not lead to the physical couplings which have non-holomorphic compo- 
nents and it is not useful for understanding the massive spectrum. In this paper we work in 
the unitary gauge and the aim of this subsection is to show how chirality can be understood 
in this gauge. 

Let us denote the affine complex coordinate on C as z and put the canonical Fubini-Study 
metric on it 

ds^ = ^"'""_ „ . (2.21) 



2 dzdz 

{l + zzf 



We now wish to turn on constant flux through C and calculate the wavefunction solutions of 
the Dirac equation. This was studied in detail in [3j whose results we use. The flux takes the 
form 

F = —^ ^dz A dz , (2.22) 

{l + zzf 

with associated gauge field 

iMz 

Turning on flux along the matter curve is known to induce chirality in the massless spectrum 
and the net chiral index is given by the integral of the flux over the curve which in this case 
is 

Xgioh.i = -^ f F = M . (2.24) 

We are particularly interested in how this chirality manifests itself in the wavefunctions. The 
massless wavefunction solutions are ^ 

l-M 

V'- = f-{z){l + zz)^ , 

1+M 

V'+ = f+{z){l + zz)^ . (2.25) 

The two wavefunction solutions are for fields of opposite charges and are the two chiralities. 
The holomorphic and anti-holomorphic functions f+{z) and f-{z) are general functions, 
which we can take to be polynomials, and therefore there are an (equal and) infinite number 
of solutions to the Dirac equation for each chirality. However chirality arises upon requiring 
that physical solutions are normalisable so that 

r r I -f 1 2 

/ dzdz^Tpl^tp± = / dzdz- t{i^M) ' (^'^^^ 

is finite. Now say M > 0, then there are no normalisable solutions for ip^ while for ip- there 
are M independent solutions corresponding to f-{z) being anti-holomorphic polynomials of 
degree up to M — 1. The opposite is true for M < while for Af = there are no normalisable 
modes. This is exactly the results of the chirality formula (I2.24p . 

Now the question we are interested in is how much of the geometry of C did we need 
to know to get these results? The chirality is determined by the convergence properties of 
zg^^± as z — )■ cxD. Say f± are taken as polynomials of order p, then the crucial object is the 
existence and position of the turning point in the function 

IP+i(l + |,|2)(-i±^)/2^ (2.27) 



1 
zg-iii± 



\z\ 



which occurs at 



FturningI — W j^,j- _ ■ (2.28) 

The existence of this turning point holds the same information as the chirality formula (j2.24p 
by counting the possible positive integer values of p. For M :§> p we have |2:turning| — ^ while 
for M = p we have l^^turningl — ^ oo. Therefore for large p we need to know the full geometry of 
C in order to count the full chirality induced by M. However for the chiral states with p < M 
we do not need to know the full geometry and for large M/p the existence of a turning point 
can be studied for values of z ^ 1 and therefore a flat space expansion is sufficient to identif; 
it. Each of the turning points can be associated to a turning point in the wavefunctions 
The extreme case being p = which is associated to a maximum of the wavefunction at the 
origin. Let us see this explicitly: in the local patch the wavefunctions ()2.25p change form and 
rather take the form of Gaussian exponentialqj 



V'- local - 


= /_(z)e 2 1^1 +... , 


V'+ local = 


= Uiz)e'T\^\' + .... 



(2.29) 

Now with the form ()2.29p we directly see which wavefunctions are normalisable or not. How- 
ever it seems that there are an infinite number of normalisable solutions since any polynomial 
of degree p would still be dominated by the exponential for large \z\. This is simply an ar- 
tifact of the fact that for large enough p the turning point where the exponential dominates 
the polynomial is forced further away from the origin and our perturbative expansion breaks 
down. Nonetheless the presence of a turning point and the associated local chirality can be 
established within the local flat patch for the wavefunctions with small p. 

2.2.2 Local Wilson lines and chirality 

In the previous section we considered an explicit global example of a matter curve. We 
showed that both chiralities solve the Dirac equation and that chirality can be understood to 
arise from the normalisability of the wavefunctions. We then showed that in the example this 
was directly related to the presence of a single peak in the wavefunctions, and that chirality 
was counted by the number of independent wavefunctions with different peak positions. 

Consider now only working in a local patch around a point in the matter curve. This is 
our setup because we must define our local patch around a particular point of interaction 
which is where the matter curve intersects other curves. The position of this point on the 
curve must be kept arbitrary in a local approach since it is only determined globally. As 
we have seen above chirality is determined by the normalisability of the wavefunctions which 
generally requires global knowledge. However we have also argued that the chirality is related 
to local turning points. We showed this explicitly for the example in section 12.2.11 and also 
showed that within this approach it is true that determining the full spectrum requires global 
knowledge but identifying a single massless mode can be done locallylfl 



®With the exception of the limiting case p — AI — 1 where the wavefunctions are a constant, though this 
exception is not present in the presence of twisting. 

The extra factor of | with respect to the form of our local wavefunctions in the exponent is because we 
have not included the twisting of the theory. 

*There may be a way to work from a local picture to a global one where identifying a wavefunction turning 
point can be associated to a global chiral mode perhaps using the approach to Morse theory introduced in [59] . 
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With these ideas in mind we can classify two types of zero modes. We fix the global 
boundary conditions and wavefunction basis and specify some point q on the curve where the 
interaction is localised. Then we define global modes as ones where the wavefunction peak 
is not located within the local patch around q, and local modes as ones whose wavefunction 
peaks within the patch. With these definitions, generally we can not ascertain information 
regarding the normalisability of the wavefunctions for global modes in the local patch and in 
particular make any statements regarding the chirality of the associated mode. Conversely, 
the definition of local modes is a chiral one: only one chirality can have a wavefunction peak 
in a given local patch. 

Let us now return to our general local E^ gauge theory. We seek a formula that will 
determine the chirality of local modes. This should depend on the local flux values as well as 
the local Higgs values. For the background Higgs and flux values ()2.2p and ()2.3p . this formula 
takes the form 

Xiocai {q") = -Re [(Mi2 + M21) mim2 + Mn {\mi\^ - |m2p)] . (2.30) 

For a given state with charges g", if xiocai > there is a localised chiral mode. In application 
to model building the formula (j2.30p presents non-trivial constraints on the possible local 
massless spectrum. The formula can be most simply derived by explicitly evaluating the 
local intersection number. The Poincae dual two-form to the matter curve can be constructed 
through the vanishing locus of the Higgs background as given in ()2.2p 



VT>[C]=^{^H)^^{^H) . (2.31) 

Then intersecting this with the flux (j2.9p gives 

F A PD [C] = xiocai VA , (2.32) 

where 774 is the volume form on the local patch. The massless and massive spectrum of local 
modes, as well as their wavefunctions, completely depends on the sign of Xiocai- We go on 
to describe this spectrum and wavefunctions in section 12.31 and in doing so will rederive the 
formula (|230]) . 

Let us discuss the relation between global and local modes, as defined above, from the 
local patch perspective. It is important to note that although for global modes the chirality 
is not determined by the local value of the Higgs and flux, once we fix the chirality, say 
from some global model, then the local profile of the wavefunction is determined and so still 
information can be gained regarding operators involving such modes. Also we should note a 
calculational advantage of local modes over global modes which is that since the bulk of the 
wavefunction is located in the patch their normalisation can be determined more accurately. 

Since the flux does not determine the particular basis of holomorphic sections it can not 
be that on a local patch we can determine whether a mode is local or global by specifying 
only the flux and Higgs, but rather it must be an additional input parameter in the local 
model. Indeed we can see this quite generally from the general forms of the local Higgs ()2.2p 
and flux (j2.3p by noting that under a coordinate transformation 

zi ^ zi + m2a , Z2 ^ Z2 — rriia , (2.33) 



Though it is not clear to the author how, or if it is possible, to formally apply these techniques to matter 
curves with flux and show generally an analogous relation between the local turning points in wavefunction 
solution and global chiral zero modes. 
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for any complex constant a, the Higgs background remains invariant. The gauge field develops 
a constant term but this is pure gauge and so can be gauged away by a phase rotation of 
the wavefunction thereby leaving its peak unchanged. Hence locally it is not possible to 
determine the position of the wavefunction peak from the equations of motion since the peak 
can be arbitrarily displaced by the parameter a. The parameter a is what we will call a local 
Wilson line, it is associated to the symmetry of the local flat metric on the patch, and if the 
metric remains flat globally so that the curve is a torus, then a is promoted to a global Wilson 
line. From a local perspective a is an input parameter that determines, amongst other things, 
whether a mode is local or global: if the peak of the wavefunction is displaced sufficiently 
away from Zi = 0, which is where the interaction point is defined, then the mode is global, 
otherwise we denote it local. We will make these statements more precise in section [3] where 
we also show that the parameter a is part of the arbitrary holomorphic prefactor present in 
local solutions to the Dirac equation on flat space, as expected from the discussion above. 

Finally, it is important to note that like other local parameters, the local Wilson line a 
can be constrained by phenomenology. A direct way to see its effect is simply to note that 
moving the wavefunction peak away from the interaction point will decrease the coefficients 
of operators involving that mode. Indeed this is the local version of the mechanism proposed 
in [10] to generate flavour hierarchies. Such a local Wilson line may indeed be useful for 
suppressing operators but as expected there is a limit to how much this can be done since the 
size of Sgut is finitelj Within a local theory a conservative bound should be that the local 
Wilson line should not displace the wavefunction peak further than the boundary of the local 
patch. We will make some quantitative statements regarding such suppression in section [3l 



2.3 The spectrum 

Having discussed the relation between local and global wavefunctions and chirality in general 
terms, in this section we describe the concrete and explicit realisation of this relation and 
how it determines the form of the massless and massive wavefunctions. 



The particular basis of the AA = 4 fields given in the explicit form for ^ (j2.1ip do not 
correspond to four-dimensional mass eigenstates and therefore we first need to go to such an 
appropriate mass eigenstates basis. To do this we follow and expand on the procedure used 
in |39 y il t [T3]. The Laplacian ()2.15p can be written as 



-AI + : 



where 



and 



(2.34) 



(2.35) 



(2.36) 



Using the D-term constraint it is simple to check that the Hermitian matrix B is also traceless. 
Let us denote the 3 normalised eigenvectors of the non-trivial block of IB as ^p with associated 









/o 







[D+,D^] [D^,Dt] 


[Dt 


\ 



^Note though that it may be possible to generate substantial suppression by deforming the local geometry, 
analogously to the findings of ^lOj where this was done with large complex-structure. 
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real eigenvalues \p which satisfy Ai + A2 + A3 = 0. Diagonalising B with a unitary matrix JI 
gives 

J^ • B • J = diag (0, Ai, A2, A3) • (2.37) 

The associated rotation of the B^ matrices is 

B+ = J^ • B+ • J* , B" = J^ • D" • J . (2.38) 

The rotated matrices D^ take the same form as (j2.1ip but with D^ — ?• D^ where 

^p" = E ^v, Dj , d; = Y. i;j d+ , (2.39) 



where ^pj denotes the jth component of the pth eigenvector. In this basis the only non-trivial 
commutators are the diagonal ones 



d;,d- 



Upq^p ■ 



(2.40) 



These commutation relations correspond to the algebra of 3 quantum harmonic oscillators. 

Let us now connect this to the general discussion of section 12.21 The massless mode 
wavefunction in the diagonal basis satisfies the equation 



D"'!' 



■^ = 



(2.41) 



where ^t = JI^^. Since we have diagonalised (j2.34p the eigenvectors ^ are simply the vectors 
with one non-vanishing entry and so there are four of them. We denote them as 



*r 








*1 







xj/. 







^, 








(2.42) 



These are the 4 massless modes that form a full A/" = 4 multiplet in the absence of a Higgs 
background and flux. Each function ipp, with P = 0, ..., 3, satisfies three equations (of which 
only 2 are independent) 

Di^o = D2'Po = D^'Po = , 
Diifi = Dt^i = D^ifi = , 
Dfip2 = D2^'2 = ^t^'i = ' 

^+(^3 = 4+VP3 = ^3-(/J3 = . 



(2.43) 



As discussed in section [2]2] chirality arises upon requiring that the wavefunctions are normal- 
isable. Consider first normalisability normal to the matter curve which is due to the Higgs 
profile. We define our local theory and patch so that the wavefunctions are strongly localised 
onto the curve and so we should see the this normalisability locally. It is simple to check, by 
turning off the flux, that ipQ and 933 are projected out by this condition leaving an AA = 2 
multiplet worth of fields. 

If we try to impose normalisability to deduce the chirality between </?i and </?2 we can only 
do this for local modes as defined in section 12.21 Global modes do not have a wavefunction 
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turning point in the local patch and so there is no sense in which they are locally nornialisable. 
For local modes the crucial objects for understanding the chirality are the eigenvalues Aj. We 
go on to discuss these but first we note a constraint on the parameters which comes from the 
local geometric picture. In this paper we require that the curves are such that their width 
is much less than their length. In terms of the background parameters, where we take as is 
expected rrii ~ Mij ~ O (1), this is the limit where i? 3> 1. The geometric regime i? ~ 0(1) is 
only consistent with the constraints on the effective theory discussed in section [2m if R± <C 1. 
Strictly speaking this is not necessarily a problem for the gauge theory as R± is a local scale 
while gravitational corrections, which are the key problem in this regime, are expected to be 
controlled by more global scales. Nonetheless it seems difficult to realise such a setup within 
a well controlled geometric regime and so we do not study this region of parameter space in 
much detail henceforth. 

In the limit R ^ 1 two of the eigenvalues, Ai and A2 scale as R and so become much 
larger than A3 which does not scale with R. Since their sum vanishes it implies Ai ~ — A2 
and we take by definition Ai < 0. The small difference between the eigenvalues Ai and A2, 
or equivalently the values of A3 is set by the background flux. We therefore can identify 3 
possibilities on their sign distribution according to the sign of the determinant of the matrix 



det B < =^ Ai < , A2 > , A3 > , 
det B > =^ Ai < , A2 > , A3 < , 
det B = =^ Ai < , A2 > , A3 = . (2.44) 

We discuss the nature of the wavefunctions for these 3 possibilities below. 

2.3.1 det B / : localised TV = 1 spectrum 

Consider the case det B < 0. In this case chirality is generated because (/?2 is projected out 
by the local criteria of normalisability. It is possible to see this by studying the wavefunction 
explicitly but there is neater approach which will also be useful for our further discussion. 
Note that in this case we see that the commutation relations ()2.40p imply that we should 
identify the raising and lowering operators of the harmonic oscillators, denoted af and a^ 
respectively, as 

a^ = iD^ , a^ = iD2 , ^3^ = iD^ , 

al = ib^ , a^ = iD:^ , a+ = W:^ . (2.45) 

With these identifications we can write the Hamiltonian as 

3 

D+D- = ^ a+ar I + diag(-Ai,0,A2-Ai, A3 -Ai). (2.46) 

1=1 

The components of ^ form 4 towers of modes and for each tower the lowering operators a~ 
annihilate the ground state and so the last term of (j2.46p gives the masses of the ground 
states. We see that this picks out the ipi function and so there is just one massless ground 
state corresponding to a single N = 1 chiral multiplet, i.e. the spectrum is chiral. The 
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particular one which is annihilated by the operators (j2.45p is 



/0\ 




pf(^ 


= 


n 


= 4^ < 


b+ip 


= 


\^) 




{Dtv 


= 



(2.47) 



These are the equations that we solve in section [3] to determine the form of the wavefunction. 
The other three eigenvectors of the diagonalised version of the Laplacian (j2.34p are of the 
same form as ()2.47p but with (/? being along the other three entries. Returning to the original 
non-diagonal basis therefore gives for the ground states 



*fi 



1 



/1\ 






('\ 



"p 



^r 



1 


^p,i 


N 


e.,2 




\lp,3/ 



V 



(2.48) 



where N is an appropriate normalisation factor. The massive excited states of the towers, 
the Landau-levels, are simply obtained by applying the raising operators which gives 



^ 



P,{n,m,l) 



iD^ 



iD-^ 



iDl 



V^;^!^(-Ai)"/2A^/2Af 
with P = 0, ...,3. The excited states have masses given b;''^'^' 



Ml 
Ml 

Ml 
Ml 



^Q,{n,Tn,l) 



l,{n,m,,l) 



^2,{n,m,l) 



3,{n,m,l) 



^p , 



M]^[-{n + l)Xi + 171X2 + 1X3] , 
M|- (-nAi + rnA2 + ^As) , 
Mk [-(.n + l)Ai + (m + 1)A2 + /A3] 
Mi [-(n + l)Ai + mX2 + {I + 1)A3] 



(2.49) 



(2.50) 



It is useful to see explicitly the mass operator in the theory since the pairing is non-trivial. 
The mass which appears in the 4-dimensional superpotential takes the form 



wM 



I I Tr 

's 



^ 



R 



^f 



R 



I I Tr 

's 



^ 



R 



T 



"^ 



R 



(2.51) 



where the indices a, b label 4-dimensional fields and here ^ are their internal wavefunction 
profiles. The mass pairing is between a representation R and its conjugate R. Note that 
by R we refer to the N = 2 partner of the state rather than the anti-particle which can 
not appear in a superpotential mass. The wavefunctions for the R states are the complex 
conjugates of the R states described above 



^ P,{n,m,l) 



^ 



R 

P,{n,m,l) 



(2.52) 



^°Note that the mass spectrum is not quite of a quantum harmonic osciUator since the ground state is 
massless. This shift in the spectrum is another way to see the requirement of the twisting of the 8-dimensional 
theory for supersymmetry. 
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but since the charge is flipped the annihilation and creation operators are interchanged so 
that now the first term on the right-hand-side of ()2.46p does not annihilate the state and so 
after performing the normal ordering there is a mass shift in the states so that 



Mi 



Ml 



M|:{-Ai,0,A2-Ai,A3-Ai} , 
M|:{-Ai,-2Ai,A3,A2} . 



(2.53) 



These are the ground states of the R modes, the excited states are created with the annihi- 
lation operators of the R states (|2.47p and have increasing masses as in ()2.50p . Note that we 
will generally drop the R and R superscripts unless explicitly required with the default state 
being R. 

The operator B~ in the mass term (|2.5ip acts as a raising/lowering operator and so 
the coupling is between states of different levels and also of different AA = 4 towers. Since 
this plays an important role in generating non-renormalisable operators in the theory after 
integrating out massive modes it is worth discussing it in more detail. The relation used to 
determine the pairing and the mass is 

^lD'-^f(„,^,0 = V^(-Al)^^f(„_l,„,;)-^/^(A2)H^(„„^_,^;) + ^//(A3)^^f(„,„,,_l), (2.54) 

i^~-^Un,m,l) = \/^^TT(-Ai)5 *^(„+i,^,,) + Vm+1(A2)5 K,(n,m+l,l) - ^Hhf' ^Mn,™,Z-l) , 



'^ 



R 



The massive modes form mass degenerate groups within which the states couple to each other. 
Let us discuss one such group as an example. First we make the external four-dimensional 
field (p and its internal wavefunction ^ manifest by using the notation 



^P,(n,m,l) — ^P,{n,m,l) (^fj) '^ P,{n,7n,l) (^0 



(2.55) 



Then using ()2.54p and (I2.5ip we find the following four-dimensional superpotential interac- 
tioi0 

W D i ($?)^ • M • $^ , (2.56) 



where 



SfQ 



^0, (0,0,0) \ 

aR 

2,(0,1,0) 
3,(0,0,1) 

v<(o,i,i)y 



M = Mk 



( 

(A2) 

(A3) 





V 



-(A2 





(As)- 



(A3)5 




(A2)5 





(A3) 
-(A2 




(2.57) 



The coupling gives the appropriate degenerate mass matrix 



M^ = Mi (A2 + As) I 



(2.58) 



which indeed matches the masses calculated using (j2.53p . Similar groupings occur for the 
other massive modes, and the pattern is replicated at each level of the tower. Explicitly for 



^^Note that here we refer to the 'physical' rather than holomorphic superpotential in the sense that the 
fields are canonically normalised. 
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some level in the tower labeled by n > we have the groupings 



^f 



SJJq 



T^3,(n+l,n,n) 
^l,(n+l,n+l,n)/ 

<t>^f ^ \ 

l,(n+l,n,n)/ 



R 



$ 



<^f 



2,(n+l,n,n) 
/,R 

0,(n,n,ra+l) 



\ 



/ 



$ 



R 



3,(n,n+l,n) 
/,R 

_2,(n,n,n+l) 



1)/ 



l,(n+l,n,n+l) 

T^2,(n,n,n) 
j,R 

'^l,(n,n,n+l)^ 



$ 



R 



\'^0,(n,n+l,n+l)/ 
\'^l,(n,n+l,n)^ 



(2.59) 



The pattern of coupling within the groups is the same as in ()2.57p : each state couples to the 
state with one more or less excitation. In integrating out massive modes therefore we inust 
make sure that they appropriately couple to each other as in the grouping above. 

The opposite case with det B > is very similar except that now the massless localised 

mode is in the conjugate representation. So if we fix the charges q"" the localised mode would 

be the one with —q'^. To summarise we showed that the sign of det B determines the chirality 

of a localised A/" = 1 mode. This then directly leads to a rederivation of the formula (|2.30p 

where we define 

det B , , 

Xiocai = -- ^ . (2.60) 



2.3.2 det B = : delocalised J\f = 2 spectrum 



In this case the spectrum changes structure. Since A3 = we see that D^ are no longer 
raising and lowering operators but rather commute with all the operators and so correspond 
to a conserved charge. The new conserved quantum number labels the Kaluza-Klein (KK) 
momentum along the curve. Now all the modes, including the massless ones, form M = 2 
representations and the spectrum is non-chiral. 

To see this we can compare the masses of the ground state (12.531) which show that as 
A3 — >• the spectrum becomes vector- like at the massless and massive levels. We now have 2 
massless modes 

D-^f = p-^R- = . (2.61) 

These are the massless modes corresponding to vanishing KK momentum. There is also a 
KK tower which generalises the equations above to 



-*R = kx"^ 



R 



^ 



R 



kK^f . 



The wavefunctions still take the form (j2.48p but now the equations for ip read 



D^ip = , DJip = 



Dfv 



-ikj^ip 



(2.62) 



(2.63) 



Now kr^ are conserved quantum numbers associated to the KK number, but note that in the 



k 



K 



local theory they are not quantised. We can define kx 
mass operator is 

D^D^if = -{kxl'^ip , 

with physical mass Mk \ kx \ ■ We can still act with the two Landau-level raising operators to 
create excited modes as before so that the modes takes the form 



(/c^) and the associated 
(2.64) 



^ 



P,{n,m,kK) 



m 



+ \n(iJJ-^m 



V^^^(-Ai)"/2 (A2) 



m/2 



* 



P,(0,0,A;if) ) 



(2.65) 
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with associated masses 

^i,„,^..^) = ^K H^ + ^ + 1)^1 + I^^H ' (2-66) 

^i,„.,..,) = < H^ + ^ + 2)Ai + IfcxH , 
^l3,,„,™,.^) = < [-("^ + ^ + 1)^1 + \^K?] . 

There is a similar grouping of massive modes as in (|2.57|) which can be deduced using (|2.62|) . 
For vanishing LL excitations, which are the only cases we will use, the pairing are 

^R ^ I '^0,(0,0,fcA-) I ^R ^ [ '^l.(0,0,fcK) ] ^R^ |'^2_,(0,0,feK) ] ^R^ I '?^3,(0,0,feA') ] 
V'^3,(0,0,fcK)/ ' V'^2,(0,0,fcK)/ ' \^l,{Ofi,kK)l ' V'?^0,(0,0,fcA')/ 

(2.67) 

The wavefunctions develop an isometry, and therefore delocalise, in this case which can 
be seen by noting that we can write 



^3^ = (^3,19., + 6,2^., j - 5 , D+= (^^3,1^.-1 + ^3,29.-,; + g , (2.68) 

where g is some function of the flux and Higgs vevs. Then using (j2.63p we have that 

(6,1^21 + 13,29^2) W?' = 'f (13,1^21 + 6,29^2 j (/? + 99 (13,1^21 + 13,2^22) f 

= \y^f [{-ikK + g) + {ikK - g)] 

= . (2.69) 

This delocalisation has of course implications for determining the normalisation of the wave- 
functions as discussed in section 13.21 Also the delocalisation shows that these modes are not 
local in any sense, indeed the presence of the massless vector pair of zero modes can not be 
determined locally and depends on the boundary conditions. 

The delocalisation also has important implications for determining the chiral nature of 
the spectrum. Although for the M = 1 modes the chirality is determined by evaluating xiocai 
at the wavefunction maximum, for the M = 2 modes, because of the isometry, we require 
that Xiocai = throughout the full matter curve and in particular also at the point of Eg. 

3 Wavefunctions and operators 

In the previous section we set up the local effective theory and showed that some information 
regarding the local massless spectrum can be obtained. Moreover the primary advantage of 
the local theory is the possibility to calculate coefficients of operators in the four-dimensional 
theory. Such coefficients come from overlap integrals of internal wavefunctions. In this 
section we calculate the wavefunctions and their overlaps as a function of general Higgs vev 
and fluxes. 

3.1 Landau-level wavefunctions 

Landau-level wavefunctions correspond to the case where Xiocai 7^ and by convention we 
choose it to be positive so that, if the massless mode is local, in the sense defined in section 
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then we know the appropriate wavefunction solves the equations (j2.47p . If the mode is 
global then the local flux does not determine the chirality of the mode, and therefore does 
not uniquely determine the wavefunction. However for concreteness we will still take the 
wavefunction to be given by (fi in (|2.42p and will include the effect of local Wilson lines. 
It may be that globally the massless wavefunction is actually ip2 in (j2.42p and this would 
correspond to the same form of wavefunction but with the charge of the mode flipped. 

In order to determine the wavefunctions for the fields as in (J2.48P we need to solve for 
the vectors ^p and for the wavefunction profile (p. The former are eigenvectors of the matrix 
IB which reads 



/o 








\ 





i (Mil + Afii) 


i (M12 + M21) 


Rfhi 





i (M21 + M12) 


i (M22 + M22) 


Rrh2 


Vo 


Rrrii 


Rm2 


/ 



Ml 



The non-trivial eigenvalues of the matrix Aj satisfy the cubic equation 

1 , , . .-. ,2 1 



(3.1) 



K + \i 



R^\mi\^ + i?V2p + ^ |Mi2 + M2if + - |Mii + Mii|^ 



+i?^Re 



(M12 + M21) mim2 + Mil (|mi 



\m2\ 







(3.2) 



Although readily obtainable, the solutions to this equation take a rather complicated form 
and so in the following we calculate in terms of the implicit solutions for Aj. A compact 
expression for the eigenvectors in terms of the eigenvalues is 



i^ 



A? 



i (M12 + M21) \i + R^m2fhi 



Mii+Mii Ai 



P2|^ |2 



,Rm2\i + \Rmi {Mu + M21] 



(3.3) 



^Rm2 (Mil + Miij 

and we define ^j = ^i/ \£,i\- There is a subtlety due to the fact that for non-generic choices 
of fluxes (for example M12 = m2 = 0) some of the vectors may degenerate in which case the 
dual expressions should be used 

/ A2 + i(Mii+Mii)Ai-fi2|^2|2 \ 

e^ = k (^21 + M12) A, + R^mim2 . (3.4) 

\RmiXi + ^Rm2 (M21 + M12) + ^Rmi (Mh + Mn)/ 

The wavefunction profile is obtained by solving (j2.47p and takes the form 

-pi\zi\'^-p2\z2\'^+P3ZlZ2+P4Z2Zl 



^ = f ( -6,2^1 + 6,1^2 



(3.5) 



where we have 



Pi 



P2 



P3 



Pa 



1 

-Mil 

2 



Rmi 



t* t* ^* ^ 

?3.2S2.3 ?2.2S; 



2,2S3,3 



1 



M22 + Rmo 



,?3,1?2,2 
'?3,1'?2,3 



?2,1?3,2, 
S2,l'?3,3 



1 



M21 + Rm2 



c^ t^ ^* t^* / 

?3,1?2,2 ?2,1?3,2/ 

/ S3,2S2,3 S2,2S3,3 



--M12 

2 



Rm-i 



,?3,1?2,2 
?3,1?2,3 



?2,1?3,2, 

(?* (?* 
?2,1S3,3 



,?3,1?2,2 S2,l?3,2, 



(3.6) 
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It is worth noting that the pure flux factors in the above expressions all cancel out in a 
charge neutral triple coupling because they appear linearly. Because of this they only affect 
the normalisation of the wavefunctions. 

Holomorphic prefactor and local Wilson lines 



As discussed in section 12.2.21 there is a local symmetry of the equations of motion which 
implies the presence of local freedom in the wavefunctions which we termed local Wilson 
lines. The values of these parameters determines the position of the wavefunction peak and 
so is the local parameter that differentiates between local and global modes as defined in 
section 12.21 We can determine how the parameter a defined in (j2.33p must appear in the 
wavefunction. Starting from (13. 5p . performing the coordinate transformation (I2.33p . and then 
performing a gauge transformation so that the background Higgs and fiux take their original 
values we find that the wavefunction transforms as 



where we define 






and 



ci 



C2 



V9 ->(/3 Cwi (a) (fwi (a) , 



g-|ap(pi|r?i2p+P2|miP+P3'^i»"2+P4m2mi) 



-pi - -Mil I "1-2 + ( -Pi + -M21 1 mi 



P2 + 2^22 j m-l + ( P3 - 2^-^12 1 "^2 



(3.7) 



(3.8) 



(3.9) 



It can be checked that the holomorphic function ip^i is a function of the particular com- 
bination of zi and 22 that appears in the general holomorphic prefactor in (j3.5p . From 
the transformation (j2.33p we know that when a ~ l/mi ~ l/m2 the wavefunction peak is 
displaced outside of the local patch and the mode should then be labeled as global. Since 
displacing the wavefunction peak along the matter curve generically is not expected to change 
normalisation of the wavefunction too much, the prefactor C^i gives an estimate of how much 
suppression of the wavefunction can be achieved by displacing it in the local patch, with an 
upper limit on a as above. Given explicit values for the background Higgs and fiux, this 
maximum possible suppression can be quantitatively analysed. 

The holomorphic prefactor is directly related to the generation structure induced by the 
flux. Indeed we can account for this by putting a generation index on the local Wilson line 

parameter ai. The actually global form of fi i — |i,2-2i + ^1,1-^2 ) will vary depending on the 
global topology of the matter curve. In the example studied in section 12.21 it took the form 
of a degree I < M polynomial, where M is the integrated value of the flux. We have some 
choice as to how we parameterise the generation dependence of the holomorphic prefactor in 
a local patch. It is possible to work directly with (pwiiai) which gives an explicit form for the 
peak of the wavefunction for each generation. 

A different basis within the generation space, is one where we take linear combinations of 
polynomials to successively eliminate all the lower powers of (—^1,2^1 + '^1,12^2) so that the 
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polynomial representative for each generation has a single power equivalent to the generation 
number 

fl {-il,2Zl + il,lZ2) = [-il,2Zl + 6,1^2) . (3.10) 

This choice of assigning the polynomial degree to the generation index was advocated in [l] 
and subsequent studies [71[8l[TTl[13] showed that it is a relevant parameterisation for some 
models of flavour such that the higher polynomials imply a smaller Yukawa coupling (and 
smaller coupling to massive modes). More generally the form (j3.10p may also have higher 
order terms, this is the case on a torus where the holomorphic basis is formed by exponentials 
[Ij that are essentially the same form as the local Wilson line factors. Whether higher order 
terms are present or not in ()3.10p can play an important role in some couplings as studied in 
section 13.41 Note that since these higher powers only give subdominant contributions to the 
normalisation integral, and so essentially decouple from the normalisation, this leads to the 
interesting possibility of suppressing couplings by suppressing the higher order powers. 

In presenting expressions we will allow for both approaches of parameterising the gener- 
ation structure. So we will consider a holomorphic prefactor of the form 

fl y-ii,2Zi + ii,iZ2J = \^-ii,2Zi + ii,iZ2J Cyji{a)(pwi{a) . (3.11) 

According to the preferred basis then different generations can be studied either by setting 
a = and using the polynomial basis, or setting the generation index / = but then 
accounting for different values for a for different generations. 

Note that, like the normalisation, the orthogonality of wavefunctions within generation 
space is a global question. One way to approach this issue is to assume global orthogonality 
by definition and thereby fix the local Wilson lines which will generically be non-vanishing 
for all the generations. The approach of taking linear combinations of wavefunctions so that 
the leading power denotes the generation number is incompatible with the assumption of 
global orthogonality in general. It may be that in some cases the two do match, such as the 
P^ example of section 12. 2. H but in general the higher order power in fi will induce mixing. 
However we can expect this mixing to be suppressed since the mixing integrals will involve 
higher powers of the Zi compared to the diagonal normalisation integrals. 

Landau-Level wavefunction normalisation 

The remaining components of ^ p as defined in ()2.48p which are required are the normalisation 
factors. Since the vectors in (12.48P are normalised the normalisation factor simply comes 
from integrating \ip\ over S. There are two components to this integral: the integration 
orthogonal to the matter curve and the integration along the matter curve. In the regime 
we are working in, where R ^ 1 the localisation of the wavefunction perpendicular to the 
matter curve is strong which means we can perform that integral as an integral up to infinity 
over flat space, this being the local approximation of 5", and expect that this does not receive 
significant corrections from the finite size of S. Regarding the integration along the matter 
curve, as discussed above, we present the result as an integral over infinity. This is expected 
to be valid for small Wilson lines and for substantial flux localisation where most of the 
integral contribution comes from the local patch. We should keep in mind that the integral 
result should not exceed the case where the wavefunction is completely delocalised along the 
matter curve, which in the homogeneous approximation is given approximately by (13.17p . 
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This latter expression may be more accurate in the case of large Wilson lines or weak flux 
localisation. Note however that, quite generally, since the localisation within the curve by 
the flux is typically rather mild the normalisation integral for a homogeneous wavefunction 
is usually not a bad approximation. Indeed the difference between the normalisation integral 
as calculated with flux localisation and without localisation was studied in appendix B of [I3] 
and was found to be relatively mild. 

In general all the wavefunction overlap integrals will be expressed in terms of a Gaussian 
integral which we define as 



/(ni,n2,n3,n4;pi,p2,P3,P4; 01,02,61,62) (3.12) 

^ni-n2 ns-m - 

/---I -Ci /Or) /On c- 
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where the evaluation is done in the homogeneous flat-space non-compact approximation for 
S. It is important to note that the result of the integral is only correct if it satisfies the 
appropriate convergence properties. 

The normalisation of the wavefunction for generation index I with prefactor (|3.1ip is then 
given by 



N? - '^ -^^m2Y- ^! 



2ki 



I - |C^,(a)i^^^^^:^ ei,2 6,1 (-ei,2ei*,ij (-ei,<2J 



2^2 / /. /v \^3/ - - \^4 



xl{ki + k3,ki +k4,k2 + ki,k2 + k3;pi +pi,P2 +P2,P3 +Pi,P4+P3; 

Ci,C2,Ci,C2) , 

where the ki are positive integers, with i = 1, ...,4, which are summed over all combinations 
satisfying the constraint X]i=i ki = I. Again we note that the result for the normalisation is 
only correct for convergent wavefunctions, while for the case of non-normalisable wavefunc- 
tions the expression simply diverges. It is also worth noting the the dependence on the Wilson 
line a in (jS.lSp drops out since the formal integral over infinity is not sensitive to coordinate 
shifts. Indeed this is why we choose to include the explicit constant factor of C^i (a) in the 
definition (fXTT]) . 

3.2 Kaluza-Klein wavefunctions 

This case is defined by the vanishing of xiocai which means that one of the eigenvalues of IB 
vanishes. We take the vanishing eigenvalue to be A3 and so Ai = — A2 and we take Ai to be 
negative. The eigenvectors of B still take the same form (j3.3p but it is simpler to give an 
explicit expression for the eigenvalues 

1 1 \ "2 

R^\mi\^ + R^\m2\'^ + -\Mi2 + M2if + -\Mn+Mnn . (3.13) 



It is informative to construct the wavefunction from the Landau- level case (j3.5p . Let us 
define (/^gauss as p.Sp in the limit A3 — )• 0. We must also take the prefactor function / = 1 
since now there is no Landau-level multiplicity and there is just one vector-like pair of states 
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at each levelo Since now Z^g commute this wavefunction is annihilated by l)j~, l)^ and 
D^. This is the wavefunction for vanishing KK momentum, while the most general solution 
should satisfy (j2.63p . We can write this as ip = (/^gaussV'KK and require 

D^ipKK = D^^KK = , DfipKK = -ik^fKK ■ (3.14) 

Solving these equation then gives the full wavefunction 

^p■. 



~-pi\zi\'^-p2\z2\'^+p-iZiZ2+P4Z2Zi+a^ zi+a^ Z2+af Zi+af Z2 /g -|^g\ 



with 
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?3,l6,2 - C3,2'^l,l C3,l'^l,2 - '^3,26,1 

ikK^2,2 -K _ ikKi2,i 



?3,1'^2,2 — ^3,2^2,1 ?3,1'^2,2 — C3,2'^2,l 



(3.16) 



Note that the i^kk part of the wavefunction is not in general a pure phase. It reduces to a 
phase if the flux vanishes which forms a subset of the cases where xiocai vanishes. Nonetheless 
the amplitude of the wavefunction always becomes flat along an isometry direction as shown 
in section 12.21 

Because of this isometry the associated normalisation integral diverges. It is still possible 
to integrate the wavefunctions in the direction orthogonal to the matter curve and what 
remains in the integral is just the volume of the matter curve. Performing this gives the 
normalisation expression 



7ri?2 Vol C 



"2 
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N^ = ^ e 4"i , (3.17) 



where we have defined 



2 



2 



ni = (pi + Pl) 6,2 + (P2 + P2) 6,1 + (P3 + P4) 6,2^3,1 + (P4 + Ps) 6,1^3,2 > 

"2 = -ih {4 + 4)+ ill {af + af) . (3.18) 

Here Vol C is the volume of the matter curve which for a near homogeneous Squt is ap- 
proximately given by Vol C ~ -kR},. This expression where the wavefunction is completely 
delocalised along the matter curve should form the upper bound on the normalisation pref- 
actor. If the localisation along the curve is sufficiently strong such that the expression (|3.13p 
gives a smaller value for the normalisation A^ then it can be trusted, otherwise, if it gives a 
larger value then the form (|3.17p is more accurate. 

It is worth displaying the wavefunction for the special case where the flux vanishes in 

which case we have 

^ = e-Ml-i-i+-2-2l' ^ (3.19) 

where we defined |7Ti| = |?Tii| + |?ti-2| . This form is important to keep in mind because we 
see that any Higgs profile localises the wavefunction onto the matter curve with a strength 
~ i?|?7i|. However, as discussed in section [2l the flux need not localise but can also repel 



""^^The absence of an arbitrary holomorphic prefactor can be seen from the fact that now we must impose 
both D^ f = and D^ f = giving two complex constraints which generically set the coefficients to vanish 
apart from the constant term. 
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the wavefunction away from the matter curve. Although it is not necessary to have the 
wavefunctions locahse to the local patch it is necessary, for the theory to satisfy its purpose, 
that the operators of the theory, or equivalently the triple overlaps, do localise. This means 
that we should require that the Higgs localisation dominates over the flux delocalisation 
including the effects of local Wilson lines. This is essentially automatic within the regime 
R^ 1 we are working mr^ 

3.3 GUT singlet wavefunctions 

The GUT singlets are by their nature more difficult to control and quantify within a local 
theory defined on a patch on Squt- These modes arise from 'bulk brane' intersections and 
so are not restricted to Scut- However the projection of the intersection curve of the bulk 
branes, on which the singlets are localised, to Squt is given by the curve on Squt as defined 
by the charges of the singlets and the associated Higgs profile {(pn)- Since the charged 
modes are localised on Squt in evaluating overlap integrals between charged and singlet 
modes we expect that replacing the bulk curve with its projection to Sgut should be a 
decent approximation. Therefore for the purpose of overlap integrals we can solve for the 
singlet wavefunctions using the general prescription outlines in the previous sections using 
the appropriate charges under the S (C/(l)^) group. 

The more serious ambiguity appears in the normalisation integral of the singlets. This 
integral should not be restricted to Squt but should go over the full bulk curve. For some 
calculations involving the GUT singlets the normalisation factor does not play a role. An 
example of this is in the relation between proton decay and exotic masses as studied in section 
14.31 However more generally the normalisation is important and forms an ambiguity which 
it seems is difficult to resolve within a local context. 

Practically we expect the normalisation of the singlets to be of the KK form ()3.17p but 
with 



Yo\C^7i[Rf''^'\ (3.20) 



2 

where i?^° ^ is associated to the global version of the local R±. Typically this is expected to 
be of order 10 for local models ^40j and so we will use ii^° ^ = 10 for numerical evaluation. 
However it is important to keep in mind that there is some uncertainty in this prescription 
for example, if the bulk brane is also local and does not stretch over the full CY volume, or if 
the flux induces chirality on the singlet curve such that the wavefunction becomes localised 
within that curve. Therefore results which depend on the normalisation may not be reliable 
to accuracy of 0(10) or so. 



^■^It is worth noting a small puzzle at this point: because of the finite width of the wavefunctions with respect 
to the Higgs curve localisation, set by the parameter R in (I3.19p . the scaling of their overlaps differs from 
the prescribed modular weights found in [55]. For example it is simple to check (see [13] for a more general 
expression including massive modes) that the physical Yukawa couplings scale as R7^R~^''^ which differs by 

the factor of 7?"^" from the usual scaling. Due to the difficulty in matching the local scales with global ones 
and the relation of R to the string coupling, it is unclear to the author whether the different scaling is a 
physical effect related to finite string coupling or due to a different parameterisation, and whether it persists 
for R±^ 3> R\\ where, as discussed in section 12.11 the local theory is unreliable. It is also worth noting that 
in [56] it was argued that the kinetic matrix could be evaluated as a residue integral over the matter curves 
rather than the full Sgut which may be related to the above issue, though the scaling of the integral given 
in [56] with R seems to be the same as would be for an integral over Sgut- 
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3.4 Wave function overlaps 

Having calculated the form of the wavefunctions we can calculate their triple overlaps which 
in turn give the associated coefficients of the operators in the four-dimensional theory. The 
relevant operator in the 8-dimensional theory is given by 

Wy= /" Tr [A A A A *] . (3.21) 



Thus, having calculated the wavefunctions for the massless and massive charged fields, cubic 
couplings in the 4-dimensional effective theory are given by integrating the triple overlap of 
the associated wavefunctions over the internal 4-cycle S. Since the wavefunctions are localised 
within the local patch so is their overlap, which means that effectively we can perform the 
integral over S as an integral over C^. 

Physical cubic couplings are determined from overlaps of normalised wavefunctions, the 
normalisation condition assuring 4-dimensional canonically normalised kinetic terms. Note 
that the canonical normalisation condition is crucial to make sense of the superpotential 
expression since by itself (j3.21|) is invariant under the full complexified gauge group [8j which 
means that it can not hold the full physical information. It is only the combined system of 
(j3.2ip and the kinetic terms that is physical. 

Wavefunctions for the charged fields can be decomposed as 

Ai = (AL ^ V'l , A2 = ./-i^ V'g , *i2 = </'iD®X^ (3.22) 

Here (/>4d are four-dimensional (super-)fields which do not depend on the coordinates on Squt 
while the internal profiles are given by tpi and x- The index / runs over all the representations 
present in the decomposition of the adjoint representation of the full enhanced gauge group 
G under the remaining gauge group after turning on the Higgs vev and fluxes. The case of 
G = Eg breaking to the MSSM gauge group is described in section [H Since we are considering 
maximal breaking to the Abelian Cartan subgroup commutant with the GUT/MSSM group, 
the generator structure of the internal wavefunctions is such that the trace in ()3.2ip simply 
leads to a selection rule stating that the charges of the three fields under the U{l)s should sum 
to zero. After accounting for this selection rule the relevant four-dimensional cubic coupling 
is given by 









/'* ..h^^^,!,:!'^ -,l,b' W.f'^^^'^'i 



1,1 I J^ll^ I 'Jlj I ^ ■,1' I J^l'^ J 

'^p,{n,m,l)Vl V2 ~ '^lp,{n,m,l)^2 ^ 

2p,{n,m,l)^ ^1 ^ ^p,{n,m,l)' 



-^tin^rnn^^'^t'+4krn,n4''^t' ' (3-23) 



We have split off from the indices {/, J, K} the generation indices {i,j, k} so that generation 
independent quantities are manifestly so. Hence a four-dimensional field is specified by fixing 
I and i. The coupling is between 3 such canonically normalised fields and we have allowed 
one of the states in the cubic coupling to be a massive state which is therefore labeled, as 
in sections [3] and 13.21 by the indices p, which denotes which of the 3 towers it is in, and n, 
m, I which label the excited state. Note that if it is a KK state the index / is replaced with 
the KK number k^- The indices i, j and k, denote the possible Landau level degeneracies 
of the states, i.e. the generation number in the case of multiple generations coming from 
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flux. Finally the factor Gijk accounts for the U{1) selection rules so that for appropriately 
normalised generators it gives 1 if the coupling is gauge invariant and vanishes if it is not. 
The triple coupling can be written as 



Y 



iI,J,K)(i,j,k) 



?p,l?l,3?l,2 "T ?p,2?l,l?l,3 "•" ?p,3?l,2?l,l 



Gijk 



I,i Tj K,k 



where the square brackets on the indices denote anti-symmetrisation. The excited wavefunc- 
tion is defined as 

'wt ) [iDo] [iD 



,/'^ = y-i; y^^) v^^3y J. 

What remains is to give the explicit expression for the overlap integral. Although the general 
expression is readily calculable we give two cases, which restrict the possible Landau-level 
excitations participating in the coupling, that can be written in a relatively compact form. 
The first case is where all the states are N = 1 Landau-levels and where one of the states 
may be a massive one with arbitrary excitation along the D^ and D^ operators. If we 
parameterise the generations by the leading power of the Zi, rather than through the local 
Wilson lines (see section 13.11 for a discussion) , there is a selection rule operating such that a 
triple coupling ^ /^'^n is generally only non- vanishing if 

i + j + k + n<m + l , (3.26) 

otherwise the integral over the wavefunctions vanishes identically. It is possible to see this by 
assigning U{1) charges to the wavefunctions where each power of Zi gives -|-1 while Zi gives 
— 1. Unless the total overlap integral has no net such charge it will vanish by the integration 
over the phases of the Zi . Note that generally, as discussed in section O the charge is assigned 
according to the leading power of Zi or Zi, but each wavefunction can also have higher order 
powers, hence the inequality sign in (|3.26p . Non-generally, if the holomorphic functions 
labeling the generations are polynomials of definite degree, as is the case when the matter 
curve is topologically P^, the inequality must be saturated. 

The selection rule implies that the leading coupling will be to states with no powers of 
the D{^ creation operator since this would add net powers of \zi\ to the integrand as it has 
positive U{1) charge. For such couplings the triple overlap is given by 

/„'-' .-^,.,-*-.t- i'^^cLcifiS, ^ \-\-\-\- f'\(i\ (''\ ("'\ M 

A"'"-"'^ " vs?.(4f' (4r''' h«hb,b,h ^''' ^"^f v*'^ w W 

e;.)'"' (-«Y.)'"" (-«&)"" (fL)" {ii,T («m)'" ('?)""■' ('?)"■"■ (*^3" ('=')" 

I {i - ki+ j - k2 + k - ks,n - ki + m- k5,ki + k2 + k3,k4 + k^; 

pi+pi+Pi,pi + pi + P2,pi +pi+pf,pi + pi+pf;ci + ci + 0^,4 + 4 + cf ,0,0) (3.27) 

where I ) denote binomial coefficients and we have introduced the expressions 



kl = ilil-pi-lMl,j+^l2(^pi-^M(2]-^lsRmi 

k'2 = ili{pi-\Mi,^+il2[-pi-lMi2)-il3Rrni, (3.28) 
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which are for i = 2, 3. Note that in (j3.27p we have taken the inequahty p.26p to be saturated. 
If it is not saturated there is an arbitrary multiplicative factor which corresponds to the 
appropriate, non-leading, term in the expansion of the wavefunctions as discussed in section 
[3j Generically this factor would be of order 1. 

In the case where the possibly heavy mode is of A/" = 2 Kaluza-Klein type the overlap 
integral takes the form 

<a'" f-«f.)'"" HQ"' (a)" ('!)"-" (*a" 



I {j - k2 + k - ks + n - k4,0, k2 + ks + k4,0;p{ + p( + pf ,pi + P2 + P2 
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where the expressions for the af^ are given in (j3.16p . We have presented the expression 
allowing only a single Landau-level creation operator to act since the more general expres- 
sion is more cumbersome as the two Landau-level creation operators do not commute. The 
analogous expression for the other Landau-level operator reads 



yn-[^2) k2=0k3=0k4=0 ^ ^ ^ ^ ^ 



■? ^2 / ^K \^ ^^ f cJ \^^ f cK \^^ n 2\^-f'4, fj2\k4 



i{j-k2 + k-k3,n- /c4, k2 + k3,k^;p{ + pi + Pi ,pi + pi + P2 ,pi + pi + pf ,pi +pi + pf; 



It is worth noting that in the case with non-vanishing KK momentum kx the selection rule 
()3.26p no longer holds for these couplings. 

Another type of cubic coupling that arises is a coupling to the conjugate representations 
R. In the case of Landau-Level modes these are, by convention, always massive (j2.53p . 
Nonetheless these couplings are important when integrating out massive modes. We will 
consider cubic couplings where the first mode, which as above may be an excited state, 
is in the R representation. As studied in section 12.3. 1|, this amounts to taking the complex 
conjugate wavefunction (|2.52p so that we are interested in, for the Landau-Levels, the integral 

^/(In.,)*^''^'"''- (3-30) 

In this case the presence of an inequality analogous to ()3.26p depends on the form of the 
holomorphic functions //. As discussed in section [3l depending on the global topology of the 
matter curve this can take the form of an exact power of the Zj, or it may be some series 
starting with a fixed power but including higher order terms. In the latter case, since the 
first term in (j3.30p includes potentially all powers of z and z, there is no selection rule at all. 
In the former case we have the selection rule 

ij^n = m + l+j + k. (3.31) 

27 



The coupling ()3.30p evaluates to 
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I ij - k2 + k - ks + m- k^ + l - kQ,i - ki+n- k4,k2 + k3 + k5 + kQ,ki + k^; 

pi+pi+pf,pi+pi+P2,pi+pi+pLpi+pi+pf-^ 

Here we have introduced 

kl = {ill) ' (-p( + ^M/i) + (ei*,2) ' (pi + ^M/2) + (ei*,3) ' Rrni , 

^2 = (li,i)' (P3 + ^^4) + (11,2)' (-pi + ^^22) + (11,3)' ^"^2 , (3.33) 

The extra coefficient a accounts for the fact that there are two possibilities of indices values 
that lead to non-vanishing results. The first is when the equality (|3.31|) is satisfied, in which 
case a = 1, and the second is where it is not satisfied and there are higher order terms in 
the holomorphic functions, a possibility that depends on the global topology of the matter 
curve, in which case a is the factor multiplying the appropriate higher order termPJ 

There are also analogous expressions for coupling to KK modes which read 

/.«»^^^-.)-^^-"^^|^t Eta) (:)(::) 
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and similarly for the other tower 
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I{j-k2 + k-k3 + n- k4,0, k2 + k3 + k4,0;pi + pf + pf ,pi + pi + p^ ,pi + pi + p^ , 
pi+pi+pf;{afY + ci + c^,{a^Y + ci + c^,{4Y^i^fy' 



^*There is a similar factor multiplying the expression (|3.27l) which we dropped because it plays a less 
important role there. 
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There are also cubic couplings that involve two massive modes and one massless one. 
These also play a role, for example in generating dimension 5 proton decay operators in the 
presence of a U{1)pq symmetry as studied in section 14.31 For the sake of brevity we will 
not display the wavefunction overlap expressions for all the combinations of massive modes. 
They can be calculated in the same way as above. Instead we present an example overlap 
which is the one used in section 14.31 This involves a cubic coupling between two conjugate 
massive KK modes (with no LL excitations) and a massless mode. The appropriate overlap 
reads 

"^■^ fc3=0 ^ ^ 

I (A; - k3,0,k3,0;pi+pi + p^,pi+pi+p^,pi+pi+pf,pi+pi+pf; 

{aff + (aff + cf , (af )' + {a^f + cf , {afY + (af )' , {aff + (aff) . (3.36) 

Since there are two massive modes involved in the coupling we should generalise the prefactor 
in the expression (j3.24p so that if we associate the AA = 4 tower index p to the field / and q 
to J we have the prefactor 

Sp,l?l,3Sq,2 ~ ?p,l?g,3Sl,2 "r ?p,2?l,l?g,3 ~ ?p,2Sij,l?l,3 T" ?p,3Sl,2Sg,l — '?p,3Sg,2Sl,l • \C>.<ii) 

4 Towards model building on the point of E'g 

The expressions for wavefunctions and their overlaps presented in the previous sections are 
quite general and hold for any gauge group. In this section we apply these results to the case 
of a point of £^8 enhancement within an SU{5) gauge group over Scut which is subsequently 
broken to SU{3) x SU{2) x U{1)y using hypercharge flux. As outlined in the introduction such 
a setup corresponds to a particularly rich theory on a local patch. The gauge theory on Squt 
has a gauge group G = Eg which is broken to Eg D SU{5)cuT x SU{5)± — )■ SU{5)gut x U{1)'^ 
by a spatially varying background Higgs vev. To see the matter content of the theory we can 
decompose the adjoint representation of £^8 under SU{5)gut x SU{5)± as 

248 -^ (24, 1) e (1, 24) (10, 5) (5, 10) (fO, 5) (5, 10) , (4.1) 

we see that we have 5 10-matter curves, 10 5-matter curves and 24 singlets. Upon breaking 
to the Cartan of 5'f7(5)_L the charges of the fields can be most concisely written in terms of 
5-component vectors (qi),- = 6ij which denote the charge under the 5 U{l)s inside f/(l)^ — 



S (?7(1)^). Then the different curves can be labeled by their charges as 



Si 



qi , 

-qi-qj, ij^j, 

(qi-qj), i^j, (4.2) 



with i = 1,...,5. To calculate the wavefunctions for a given matter curve we contract the 
associated charge vectors with the background Higgs and fiux ()2.2p and (j2.3p which gives the 
effective Higgs and fiux felt by the modes on the matter curve as in (I2.14p . The resulting 
wavefunctions overlaps are then calculated as in section 13.41 which give the coefficient of the 
appropriate cubic coupling in the four-dimensional theory. 
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m, , ..... ...,^. 



Therefore in specifying the local theory the input parameters are the background Higgs 
flux M?- along the S (C/(l)^) x C/(l)y generators, the geometric quantities i?|| and R±, 
and the local Wilson line values ci and C2 for each wavefunction. Once these are specified the 
local form of the wavefunctions is determined and the resulting overlaps can be calculated. 
As well as determining the overlaps, as discussed in section [2?2t in the absence of local Wilson 
lines the local values of the fiuxes are connected to the chirality of the localised modes. 

The variety of possible local flux and Higgs values, and the intricate relations between 
these quantities and the spectrum and operators of the effective four-dimensional theory leads 
to a rich and complicated phenomenology covering a wide range of phenomenological issues. 
A comprehensive study of the various possibilities and their phenomenological implications is 
beyond the scope of this paper. Instead, at this stage, the aim is to introduce the important 
concepts and provide the necessary tools to study these theories. To this end in this section we 
will study an example model which is not chosen for particularly attractive phenomenological 
properties but rather simply to allow for explicit, rather than abstract, discussions of the type 
of analyses that can be performed for these models. 

Some tools for analysing such models are provided at: 



http : //www . cpht . polytechnique . f r/cpht/palti/pointES . html 



The web page takes the background values for the flux and Higgs vevs as input values. 
It applies (j2.30p to the different effective values for all the representations calculating the 
resulting chiral spectrum in the local patch ignoring the possible local Wilson lines, i.e. it 
provides the chirality induced in the case where the Wilson lines vanish so that we can be 
sure that the local patch captures this correctly. It also calculates the effective Higgs and 
flux felt by each representation according to its charges and there is a link to a Mathematica 
worksheet which calculates the resulting wavefunctions and overlaps for 3 such matter curves. 

A condition which we impose on the Higgs field vev is that locally the gauge rank does 
not enhance beyond 1 rank except at the point of Eg at the origin. This is not a necessary 
condition for general models, for example it is possible to turn off matter curves all together, 
or to enhance by rank 2 over a matter curve. However for simplicity we would like to require 
the above property, which also ensures that the triple overlap of any wavefunctions, in the 
presence of just a Higgs background, is localised at the origin. There are of course many 
values of Higgs fields that satisfy this property but a simple way to ensure it is to take the 
following ansatz 

m? = e^^^"/^ , m^ G Rationals / . (4.3) 

Since the fifth roots of unity form a cyclotomic extension of rank 4 over the rationals it is not 
possible to set equal any two curves for rational non-vanishing values of the 1712 (other than 
using the single linear relation which corresponds to the tracelessness constraint). Of course 
in practice there may be approximate cancellations which could lead to near vanishing Higgs 
profiles along some directions we should be aware of. 

The chosen example model is defined by the local values for the Higgs and Flux back- 
grounds given in table [H With the input values specified we go on to discuss aspects of this 
model as explicit cases of more general properties of these type of local models. 

Before proceeding it is worth fixing the local geometric parameters R± and Rn within the 
local theory. Given that the Higgs vevs taken are O (1), we choose to study the geometric 
parameter space by fixing R± = 3/2 and Rn = 3 which is consistent with the constraints 
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Higgs 


Value 


Flux 


Value 


m\ 


1 


m\ 


4 
5 


2 

nil 


e^TTj/s 


7712 


-1 


ml 


g47rj/5 


Mh 


11 

5 


m\ 


„67ri/5 


Ml 


14 
5 


m\ 


g87ri/5 


Mf, 


_11 
5 


7TT-2 


3 

5 


Mh 


8 
5 


7TT-2 


Y 
5 


Mil 


y 

5 


m^ 


1 







Table 1: Table showing the non- vanishing Higgs and flux values for the example model. The 
super-indices denote the S (C/(l)^) x U (l)y generators. Note that the D-terms are implicitly 



solved so that M22 



-Mft. 



Representation 


Xlocal 


f(3, 2)1/6 (3, 1)_2/3©(1, 1)1) 0qi 


+ 1 


((3,2)i/6©(3,l)_2/3©(l,l)i)®q2 


+ 1 


((3,2)i/6©(3,l)_2/3©(l,l)i)0q3 





((3,l)i/3e(l,2)_i/2)c5(-qi-q3) 


-1 


((3,l)i/3©(l,2)_i/2)®(-q4-q5) 


-1 


(3,l)_i/3(8)(-qi-q2) 





(l,2)+i/2®(-qi-q2) 


+1 


(3,l)_i/3(g)(-q2-q3) 





(l,2)+i/2®(-q2-q3) 


-1 



Table 2: Table showing the chirality for local modes. 

outlined in section [2j Note that, as emphasised in section [2] and in |J^, the local value for 
i?_L would have to differ by a factor of O (10) from the global value which sets the Planck 
scale; this being a universal issue with local gauge theory models. 



4.1 The massless spectrum 

The first property of the local models to study is the massless spectrum. Since the model is 
defined on a local patch within Squt it can not provide full information regarding this but 
rather is connected to it through the discussion presented in section [2^21 We choose out local 
Wilson lines so that the local modes include those given in table [21 By definition these modes 
have vanishing local Wilson lines implying their wavefunctions peak at the origin and their 
associated local chirality is determined by the background Higgs and fluxes as presented in 
table [1] applied to the formula ()2.30p . The local chirality for these modes is shown in table [2j 

The local spectrum is compatible with exactly the matter content of the MSSM, with 2 
generations living on one of the 5 matter-curve and similarly 2 generations on one of the 10 
curves. All the other modes may have local Wilson lines and so their chirality is unrelated 
to the local values of the Higgs and fluxes. However, their local values for Xlocal are non- 
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Field Name 


Representation 


Multiplicity 


(Qiec/ieSi) 


((3,2),/6e(3,l)_2/3e(l,l)i)»qi 


1 


{Q2 e c/2 e E2) 


((3,2),/6e(3,l)„2/3e(l,l)i)®q2 


2 


PieLi) 


((3,l)i/3e(l,2)_i/2)®(q4 + q5) 


1 


(1)2,3 e ^2,3) 


((3,l)i/3®(l,2)_i/2)®(qi+q3) 


2 


Hu 


(l,2)+i/2«>(-qi-q2) 


1 


Hd 


(I,2)_i/2®(q2 + q3) 


1 


El 


(i,2)_i/2»(qi + q5) 


1 


E2 


(1, 2)^1/2 ^(-qs-qs) 


1 


- 


((3,2)1/6 02 X (1,1)1) 0q2 


1 


- 


((3,2)_i/6©2x(l,l)_i)0(-q5) 


1 


- 


(3,l)i/3 0(q4 + q5) 


1 


- 


(3,l)_i/3®(-q2-q4) 


1 


Xi 


(I,l)o0(-qi+q3) 


1 


X2 


(l,l)o(»(-qi+q2) 


1 



Table 3: Table showing the massless spectrum of the example model as determined from a 
semi-local approach. 



vanishing which means that although there need not be necessarily any net chirality over the 
matter curve, or any massless modes at all, we also know that many modes can not be of KK 
type since that requires xiocai = over the full matter curve. 

It is important to note that a motivation for considering models where most of the modes 
are local is that the presence of their wavefunction peak at the interaction point maximises 
they size of the coefficient of the operators involving them. This is particularly important in 
the context of recovering sizable Yukawa couplings. 

We would like to embed the chirality of the local spectrum into a global model. We can 
realise this, at least within the semi-local approach, by considering the semi-local spectrum 
presented in table [3l The presence of additional non-MSSM modes is required by a consistent 
embedding of the spectrum into Squt- The constraints of this embedding are calculated in 
the appendix where also the explicit embedding of the spectrum of table [3] into Squt is 
presented. We have also added a massless GUT singlet field X^j The spectrum of the 
example model is in many ways typical of the type of models studied in semi-local F-theory 
GUTs, for example in [411143] . It has additional U{1) symmetries with selection rules, in 
particular this model exhibits a U{1)pq symmetry which forbids dimension 5 proton decay. 
It has exotics which are associated to the presence of such a symmetry which can only be 
lifted by giving a vev to the singlet X|^ An important aspect by which it differs from more 



^^The net chirality of the GUT singlets is not constrained in the same way as the charged matter since it 
is sensitive to the pull back of the flux to divisor intersections away from Sgut and therefore require a full 
global model to determine. Some local constraints of the chirality of the singlets that arises just from the 
component of the flux which pulls back to curves on Sgut was studied in [411142] . 

^''There are other exotics in table[3]that are unrelated to the presence of a U(1)pq which can also be lifted 
by giving appropriate singlets a vev. 
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advanced F-theory models is the rank of the up-type Yukawa matrix being 2 rather than 1 
because to achieve a rank 1 matrix requires local monodromies [5tl35j. 

Having defined the massless spectrum of the local theory we can calculate the associated 
operators. Over the upcoming sections we do so for various type of operators that are of 
phenomenological interest. It is important to emphasise that primarily due to an ambiguity 
in the normalisation integral of the wavefunctions, as discussed in section [3l the numbers 
quoted are only expected to hold to within factors of 0(1) (or perhaps 0(10)). 

4.2 Yukawa couplings 

In this section we study various aspects of Yukawa couplings. We begin by considering the 
Yukawa coupling to Hu, for example YhuQiU2- ^11 of the fields participating in this coupling 
are localised within the patch and so this coupling serves as interesting case for studying the 
dependence of the magnitude of couplings on the approximations used in the local effective 
theory. Using the formula for the integral ()3.27p gives the coupling 

YhuQ.U, ~ 3 X 10-3 , (4.4) 

where here, and henceforth, we quote the absolute value of the coupling. Our conventions for 
the normalisation of the wavefunctions is that, for a given wavefunction, out of the Landau- 
Level normalisation (|3.13p and the Kaluza-Klein normalisation (|3.17p we use whichever one 
gives the smallest value for the normalisation factor A^. The motivation being that if the flux 
localises the wavefunctions significantly within matter curve then the LL normalisation will 
be more accurate than the KK one which is approximate. 

It is clear that the magnitude of the coupling (14. 4p is far too small to reproduce the top 
quark Yukawa coupling which must be close to 1 at the GUT scale. Since the kinetic terms 
are canonically normalised there is no overall scale that can be used to reconcile the two 
values. This problem is quite general and is associated to the difference in localisation of the 
wavefunctions perpendicular and parallel to the curves which means that the wavefunctions 
normalisation integral is too large compared with the overlap integral. Generally couplings of 
O (1) are associated with the large flux limit, Mij ^ 1, which implies that the wavefunctions 
become strongly localised withing the curve, or the geometric limit Rn <C 1. Neither of these 
limits are compatible with the validity of the local theory, however it is plausible that for some 
regions in the valid local parameters space that are closer to this limit a more substantial top 
Yukawa coupling can be induced. 

The down-type Yukawa coupling is calculated to be 

YhoQiD, = 3 X 10-3 . (4.5) 

Note that an important aspect of studying the point of Eg is that, as we have shown, both 
the up-type and down-type Yukawa couplings can be calculated within the same local theory 
allowing for a determination of their relation and consequently parameters such as tan/3. 

Non-normalisable Yuka^va couplings 

The generation of Yukawas for the lighter generations is less clearly understood. The reason 
being that the presence of a normalisable Yukawa coupling for one generation, associated with 
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at least a rank 2 enhancement of the gauge group over a point, does not imply non-vanishing 
Yukawas for the lighter generations [8','9','36] . Therefore some additional features are required 
to generate the other Yukawas and these depend on how the generations are realised. In the 
case where the generations are Landau-level multiplicities coming from a single matter curve, 
the most natural possibility is that there are a number of intersection points along the curve 
so that at each one a Yukawa coupling is induced. This possibility is the most natural from a 
top-down perspective, however it is less natural when we account for the observed hierarchies 
of quark and lepton masses in the SM. It is not clear how and why this should naturally 
arise within such a construction. Nonetheless it is a viable possibility and we refer to [10] for 
further studies of this scenarioPl 



Another possibility, studied in [8l l36l[TT] . is that non-perturbative effects induce a non- 
commutative deformation of the local theory which does induce a rank 3 Yukawa coupling 
associated to a single enhancement point. This itself has difficulties particularly in embedding 
into a global realisation, and in generating Yukawas that are not too small for the lighter 
generations given that they have a non-perturbative origin. 

A qualitatively different approach is to associate each of the generations to a different 
matter curve. This way it is possible to control the structure of Yukawa couplings using the 
additional U{1) symmetries naturally present in breaking E^ to SU{5) [46-48J. In this case 
some of the Yukawas arise as non-renormalisable operators in the four-dimensional effective 
theory after integrating out heavy modes. Such a setup is particularly interesting to study 
within the local context since all the important couplings can be calculated. For the rest of 
this section we will study this in detail since it is also relevant for generating an approximate 
rank 1 up-type Yukawa coupling, as opposed to the rank 2 in the example model. l£| 

In the example model of table [3] a non-renormalisable Yukawa coupling which is gauge 
invariant is given by 

Y^^X2HuQiUi . (4.6) 

This is an interesting case as it amounts to a four-dimensional version of the monodromy 
structure typically used in F-theory GUTs to generate a rank 1 up-type Yukawa matrix. The 
coupling is induced by integrating out massive modes from the renormalisable couplings 

^X2QiQ2^2<5lQ2 +YhuQ-2UiHuQ2 Ui , (4.7) 

where the superscript LL denotes massive Landau-Level states as studied in section ()2.3.ip . 
Because of the fact that the mass operator of massive LL also changes the level and so the 
wavefunction these couplings are non-trivial in the sense that they require an understanding 
of the holomorphic part of the wavefunctions beyond the leading expressions (I3.10p . This is 
because, as seen in the model of table [3l there is only one generation on each of the curves 
carrying Hjj, Qi and Ui which means that the leading parts of the holomorphic prefactor 



^^Note that such a scenario could be studied within a local patch though the distance between the intersection 
points would push the viability of the local theory to the edge. 

^*It is important to note that as we have seen with the top Yukawa the degree and form of localisation 
plays an important role in the magnitude of the coupling which naturally raises the possibility that the 
Yukawa structure could arise from the localisation properties of the generations. This has also been suggested 
in a large number of other works, see |49II5U| for old and recent references. Indeed we have seen that an 
accurate determination of the top Yukawa within a local theory hints at point-like localisation for some of the 
wavefunctions involved. Either way the localisation properties of the wavefunctions clearly play an important 
role in determining the magnitude of the operators and it certainly deserves further study within local theories 
where the associated parameters can be more accurately determined. 
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of the wavefunctions is constant. However this part can not contribute to the couphng ()4.6p 
because in (j4.7p either Q2 ^ or Q2^ must have a non-trivial dependence on net powers of Zi 
which would mean that the appropriate cubic wavefunction integral would vanish. 

Therefore, in evaluating a non- vanishing contribution to ()4.6p from integrating out massive 
LL, there is an ambiguity relating to the coefficient of the subleading holomorphic power in 
the wavefunction. As discussed in section [3l there are two unknown factors in determining 
this coefficient, the first being if there are subleading powers present at all which depends on 
the topology of the matter curve. For the example studied in section 12.21 of P^ there were 
no such powers as they would lead to non-normalisable wavefunctions, while in the case of 
a torus such higher order terms are present. Assuming that such terms are present there is 
then an arbitrary coefficient multiplying them which only gives a sub-leading contribution to 
the normalisation integral and so is relatively unconstrained though expected to be of order 
one. We will calculate the contribution to (j4.6p that arises from a linear holomorphic term 
in the wavefunction prefactor of X2. This is expected to be present since the wavefunction 
for X2 is the pullback of a bulk wavefunction to the GUT brane anyway and so is expected 
to have higher order powers, further since we have little control of the normalisation of the 
singlet wavefunction and have to estimate it the ambiguity in the coefficient of the linear 
holomorphic term is absorbed into the normalisation. Therefore, we should keep in mind 
that the result can only be trusted to within factors of O (10) or so. 

From the grouping of LL calculated in (|2.59p we see that the appropriate heavy modes to 
integrate out are the pairings 

{Q2 Ml }= 1*2,(0,0,0)' ^3,(1,0,0)/ '1^3,(0,0,0)' *2,(1,0,0)| ' (^•^) 

These are expected to give the leading contribution to the non-normalisable operator. La- 
beling the coupling to the first set of heavy modes with superscript 1 and to the second with 
superscript 2, taking the coefficient of the linear term in the singlet wavefunction to be 1 and 
using the approximate singlet normalisation discussed in section 13.31 the formulas presented 
in section 13.41 give 

^kQ.Q. ^ 7 X 10-^ , n,Q,v. ^ 5 X 10-^ , 

^iQ.Q. - 10"' ' "^kQ.V. - 10-^ . (4.9) 

The contribution from the heavy modes ()4.8p to the coupling ()4.6p is therefore calculated to 
be 

^^^ ~ MK^f^l V^^Q^Q2^^uQ2Ui + ^X2QiQ2^HuQ2Ui) j^^ ■ (4-10) 

Even though this figure may not be accurate to within factors of 0(10), it is still interesting 
to note that it is suppressed with respect to the tree-level coupling (j4.4p . This may have 
implications for example for the models presented in [46j which rely on the two being equal 
and induce additional suppression of the non-renormalisable coupling through a small singlet 
vev. Of course it is plausible that for some range of geometric parameters this may hold but 
for this particular model, at the level of accuracy to which we can work, it seems that there 
is additional suppression coming from integrating out the heavy modes. 

It is important to keep in mind though that there is a whole tower of massive modes 
which must be integrated out, whereas we have calculated just the leading contribution. Since 
higher mass modes also have higher powers of Zi in their wavefunctions their contributions 
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are suppressed by both the higher mass and the smaller coupling to the massless modes. We 
therefore expect the leading contribution to be a decent approximation. 

Hypercharge flux and b-T uniflcation 

An important piece of evidence towards grand unification is that at the GUT scale the 
masses of the bottom quark and tau lepton unify (see [51] for a modern analysis). There is 
also approximate unification for the lighter generations. Given that hypercharge flux modifies 
the wavefunctions for the bottom quark and tau lepton it is important to consider whether 
these mass relations are maintained. Indeed it was suggested that the slight non-unification 
for the lighter generations may be sourced by the hypercharge flux [llj. Here we are more 
interested in whether the accurate unification for the heaviest generation is maintained and 
so calculating the appropriate coupling we find 

YHnQiD, = 3 X 10-3 , Yh^e.l, = 4 X 10-3 . (4.11) 

To within the accuracy we are working with this is a decent match which shows that it is 
possible for the two to match. However it is important to state that the two are not the 
same and generally in the presence of hypercharge flux the two need not match which places 
a constraint on the magnitude of hypercharge flux near the bottom Yukawa point. 

4.3 Exotics and Dimension five proton decay 

An important constraint on GUT theories comes from proton decay operators. In this section 
we will study dimension 5 proton decay operators that are induced by integrating out mas- 
sive modescj The spectrum was chosen so that it exhibits a U{1) symmetry which forbids 
dimension 5 proton decay operators. This is a common idea in F-theory GUT model build- 
ing but such a symmetry, combined with the use of hypercharge for doublet-triplet splitting, 
implies that there are exotics in the massless spectrum whose mass is protected by the same 
symmetry |45 p 41 t l43 p52H54j . In this example these exotics are marked as Ei and E2 in table 
[3l It is possible to lift these exotics by giving a vev to the singlet Xi but this will then induce 
dimension 5 proton decay operators. There is therefore a natural tension between raising the 
mass of the exotics, for example to comply with gauge coupling unification, and suppressing 
dimension 5 proton decay [MKS^. Studying this generic situation requires understanding 
quantitatively the coefficients of the exotics-singlet coupling and of the dimension 5 proton 
decay operator. An initial study in understanding the generation dependence of one of the 
couplings that features in the calculation, the Yukawa-like coupling of the massive Higgs 
triplet, was performed in [13j. In this section we will extended this to a full calculation of the 
dimension 5 operator. Although of course the results are not general, and can be modified, 
for example by the presence of monodromies around the up-type Yukawa, this will serve as 
an interesting case study of such a process. Further if we are only interested in the ratio of 
the coefficient of the dimension 5 operator to the coefficient of the cubic coupling used to 
lift the exotics the normalisation of the singlet matter curve cancels thereby increasing the 
accuracy of the calculation. 



Note that the spectrum in table [3] admits dimension 4 proton decay operators though this does not affect 
the use of it as a case study for dimension 5 operators. Although the model is not proposed as phenomeno- 
logically viable in itself, for this and other reasons, it is worth noting that assigning all the generations to the 
curve q4 + qs will forbid these operators. 
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The dimension-five proton decay operator arises after integrating out massive states in 
tlie cubic interactions 

W D YuT^''QiQ2 + YdT^'^QiLi + YxXif^^'f^'' , (4.12) 

which, after integrating out the first massive modes, lead to the operator 

^ ^ ^^Zl^^i^i^^Qiii . (4.13) 

Here T^^ denote the massive KK-type modes, following the spectrum presented in table 
[21 of the Higgs triplets. Note that when integrating out the modes there are 4 possibilities 
according to the pairings (I2.67p . corresponding to choosing the KK mode to be in the first or 
second tower for the Higgs up and Higgs down fields. All of the possibilities contribute to the 
higher dimensional operator but for brevity we will initially consider only the contribution 
coming from the choice where both T^^ and T^^ are in the 2nd tower (and therefore T^^ 
and T^^ are in the first tower. We also have the coupling of the singlet to the exotics 

W D Yx,E,E,XiEiE2 . (4.14) 

We can calculate all these couplings within the local theory using the framework presented 
in section 13.41 It is simple to check that for the exotic Ei the local chirality is opposite to 
the global chirality which means that for consistency we have to turn on substantial local 
Wilson lines in its wavefunction which we do by taking a = 1 for that mode. Taking the KK 
momentum in ()3.16p as kx = 1 gives the couplings 

y„ ~ 5 X 10"^ , Fd ~ 2 X 10""^ , 
Yx ^ 2,x 10-^ , Yx^E^E2 ^ 5 X 10^3 ^ (415) 



which implies 

YuYdYx 



10"*^ , (4.16) 



Yx^E^Ei 

while we also have 

YhoQ.d,Yh^q,u, ^ 10-^ . (4.17) 

An interesting thing about (j4.16p is that the proton decay operator is substantially suppressed 
compared to the operator that lifts the exotics, even further, by a factor of 10~^, than would 
be expected given the suppressed Yukawa coupling values ()4.17|] . Although this is a toy model 
it is encouraging to note that substantial further suppression factors are possible and could 
possibly alleviate the tension between proton decay and gauge coupling unification discussed 
above by allowing for large exotics masses. 

5 Summary 

In this paper we studied operator coefficients in local theories of F-theory GUTs by calculating 
wavefunctions and their overlaps. We solved for wavefunction profiles and overlaps in general 
linear Higgs and flux backgrounds for general gauge groups. We then applied these results to 
a model based around a point of E^ enhancement and calculated coefficients of some of the 
operators in the theory. 
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The model for which we calculated the coefficients is just a toy model used to illustrate 
the type of calculations that can be performed. The aim of this paper is to provide some of 
the tools and framework that allow for studying more generally the type of theories that arise 
near a point of Eg. The interplay between the local values of the metric, flux and Higgs field, 
and the local spectrum of operator coefficients is far too rich to have been systematically 
studied. For example, although we developed the formalism, we did not use local Wilson 
lines in the example model to suppress operators. Even the type of operators we studied are 
only a small set of those based at a point of Eg, and other operators can be studied using the 
methods and results of this paper. It would be very interesting to study more methodically 
the possible phenomenology that can arise from these models in future work. 

However, even within the current study, there are some general phenomenological aspects 
worth mentioning. It seems that it may be a challenge to recreate an 0(1) top Yukawa 
coupling within a valid local effective theory. Indeed in the example model we found that 
it is too small. Although this does not have to be so always, generically the reason is that 
the coupling is suppressed by the size of the local patch, effectively by its modular weight 
under the Squt modulus. We suggested that this may point towards point-like, rather than 
curve-like, localisation of the fields involved which can be achieved by strong fiux amongst 
other things o A more general phenomenological conclusion is that quite generally factors 
that naively may be of 0{\) can easily be a few orders of magnitude smaller. A particular 
examples where this turned out to be useful is in calculating dimension-five proton decay in 
the presence of a (broken) U{\)pq where we found that the non-renormalisable operator was 
more suppressed than expected. 

A primary theme of the work presented is a study of the constraints on the local theory 
which turn out to be non-trivial. For example we argued that generically the local approach 
of using an 8-dimensional Higgsed theory to describe an intersecting brane setup, is not valid 
over the full Sgut if the normal directions to it in the CY are large, the so called decoupling 
limit. Indeed this was taken as a practical motivation for working in a patch within Sgut- We 
also studied in detail the relation between the local values of the flux and the local massless 
spectrum which imposes restrictions on both. 

There are a number of phenomenologically motivated improvements that can be made to 
the simplest models studied in this paper. An important one is the incorporation of local 
monodromies, perhaps using the tools presented in |35j . This is likely to require numerical 
methods in solving for the wavefunctions since complete analytic solutions to the type of 
differential equations that arise in this case are not known. 

Perhaps more radically it may be that the size of operator coefficients, or natural global 
embeddings, require that we go to regions of parameter space that are beyond the validity 
of the local theory itself. In particular this applies to the case where the component of the 
metric of the normal direction to the brane is larger than the parallel ones. In that case a 
more appropriate theory would seem to be a 10-dimensional one with a product gauge group, 
ie. separating the bulk branes out rather than treating them as small deformations of the 
GUT brane. It is not clear how much it would be possible to calculate in such a theory for 
a general bulk geometry and whether the local approach can still be applied. On a similar 
theme it would be interesting to see if there is some sense of a local theory that can be applied 
to Heterotic constructions where the localisation of the wavefunctions would be done purely 



^'^In [35] a proposal for localising point-like matter was made using local monodromies, it would be interesting 
to see if such localisation could be used to enhance the absolute value of the top quark coupling. 
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by flux effects. In particular if it is possible to localise using fluxes onto the analogue of a 
point of £"8 perhaps using the line-bundle constructions of 



As discussed in the introduction the local approach to model building has a number of 
advantages, in particular to do with calculability and generality within the landscape. In this 
work we have taken some steps towards this aim for a simple class of models. There is much 
work to be done if we are able to reach the full potential of a phenomenologically viable local 
model which recreates, at least approximately, the operator coefficients that we observe and 
constrain in experiments. And of course it may be that the real world does not correspond to 
a local model at all. But it seems to be a direction worth pursuing in that it has the potential 
of forming a relatively direct route to quantitative comparison of a large and general class of 
models with a large experimental dataset. 
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A Semi-local constraints on matter spectrum 

In this appendix we calculate some constraints on the restriction of hypercharge flux coming 
from the embedding of the matter curves into the full Squt- Given a spectrum of fields 
within a local patch on Squt as discussed in section 12.21 the constraints on the full Squt 
derived in this section can be used to understand what the spectrum away from the local 
patch must be like in order to be compatible. 

The techniques used to derive the homology classes of, and therefore the possible flux 
restrictions to, the matter curves are described in detail in |43) . to which we refer for further 
details and here we will simply state the results. Following precisely the methodology and 
notation of j43j we start from a fully split spectral cover 



Cio = (aiV + 06^7) ia2V + ajU) (a^V + agU) (a^V + agU) (a^V + awU) = . (A.l) 

Here the aj are some as yet undetermined coefficients that are functions on S. The spectral 
cover can also be written as 

Cio = boU^ + b2V^U^ + b-iV'-^U^ + b^V^U + 65^^ = , (A.2) 

where we know that the 6j are zero sections of the bundle r] — ici |44lH5] . Here ci is the first 
Chern class of the tangent bundle of Squt and r] = 6ci — t with —t being the first Chern 
class of the normal bundle to Squt- Comparing ()A.2p and (jA.ip then implies that the a/ 
are sections of bundles as shown in table HI Here X{7,8,9,io} are unspecified and we define 

X = X7 + X8 + X9 + XlO- 

To obtain the homology classes of the matter curves we need to solve the constraint 61 = 
which we do by setting 

ai = — (020809010 + 030709010 + 04O7O8O10 + 05O7O8O9) , 

06 = 07O8O9O10 , (A.3) 
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Section 


ci (Bundle) 


ai 


?? - ci - X 


02 


-Cl + X7 


03 


-Cl + X8 


04 


-Cl + X9 


05 


-Cl + Xio 


ae 


v-x 


a? 


X7 


as 


X8 


«9 


X9 


OlO 


Xio 



Table 4: Table showing the first Chern classes of the line bundles that the aj are sections of. 
The forms X{7,8,9,io} ^^^ unspecified and we define X = X7 + Xs + X9 + Xio- 



which also imposes the homology constraint 

f/ = 2x • 



(A.4) 



With this choice it can be checked that there are no exotic non-Kodaira singularities on Squt 
over any of the curves Oj o 

The matter curves are given by the equations 



P5 



hlhA - h2hh + bobj = 



(A.5) 
(A.6) 



Which, after imposing ()A.3p . give the curve homology classes as specified in table [5l The 
table also shows the possible flux restrictions over Sgut for the matter curves which in 
turn constrain the possible chiral spectrum. Given the flux values Ni and Mj the resulting 
spectrum is determined as 



for the 5-matter curves and 



"(3,l)-i/3 - 


" "(3,l) + l/3 


= Ms, 


"(l,2)+i/2 - 


- ^(l,2)_l/2 


= Ms + iV, 


^{3,2)+i/e - 


^{3,2)_i/6 = 


= Mio , 


'^{3,l)-2/3 " 


'^(3,l)+2/3 = 


= Mio-iV, 


^(l,l)+i 


-^(1,1)-1 = 


= Mio + iV, 



(A.7) 



(A.8) 

for the 10-matter curves j^i 

The example model of section H] is realised explicitly by taking the non-vanishing flux 
values presented in table [H 

^^Note that the relations (|A.3[I solve the constraint 61 = also with an additional arbitrary section c 
multiplying each of the right-hand sides. However if c has a zero section then it induces a non-Kodaira type 



singularity on Sgut and so we should set c = 1 which implies the form (|A.4[l . 

^^It is worth noting when comparing the local and global values of the fluxes that the integrated flux 
values here M and A'' are combinations of the hypercharge and U{1) fluxes such that M5 = M'^^^' — iN^ , 



Ml 



M^(i) + iAT^, TV : 



iN' 
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Matter Charge 



Equation 



Homology 



N 



Mu(i) 



5i 


qi 


q2 


52 


qi 


q3 


53 


qi - 


q4 


54 


qi - 


^5 


55 


q2 


q3 


56 


q2 


q4 


57 


q2 


qs 


58 


q3- 


q4 


59 


q3 


qs 


5 10 


-q4- 


qs 


lOi 


qi 




IO2 


q2 




IO3 


qa 




IO4 


q4 




IO5 


qs 





asogaio + aia^aif) + a^asag 

o-20-90'io + dAdjaiQ + a5a7a9 

a2a8aio + a^ajaio + 050703 

020309 + 03O7O9 + 04O7O8 

O2O8 + O3O7 



-ci + X8 - 

-C\ + X7 ^ 
-ci + X7 ^ 
-Cl + X7 - 



X9- 
X9- 

X8- 
-X8 



Xio 
Xio 
Xio 

-X9 



0209 - 
O2O10 
03 09 ' 
O3O10 
04010 



0407 

-0507 

0403 
-0508 
- 0509 



-Cl + X7 + X8 
-Cl + X7 + X9 

-Cl +X7 + X10 

-Cl + X8 + X9 

-Cl -H X8 + XlO 

-Cl -H X9 + Xio 



-N^ 

~Ng 

-iVio 

A^7 + ^8 



Nr- 
Ns- 

Ng- 



-Ng 

iVio 

'Ng 

NlQ 



01 

02 

as 

04 

05 



-Cl 

-Cl - 
-Cl - 
-Cl - 

-ci4 



HX 

-X7 
-X8 
-X9 

XlO 





iV7 

Ns 

Nq 

Nio 



Ms, 
Ms, 
Ms, 
Ms, 
Ms5 
Ms, 
Ms, 
Ms. 
Ms, 
Msio 
Mioi 
Mio, 
M103 
Mio, 
M105 



Table 5: Table showing curves and flux restrictions on Squt in the absence of monodromies. 
The fluxes are constrained so that N-j + N^ + Nq + Niq = 0. 



Flux 


Value 


M5e 


+1 


M5, 


-2 


^5io 


-2 


Mioi 


+1 


Mio, 


+3 


^105 


-1 


iVr 


-1 


iVio 


+1 



Table 6: Table showing the flux values used to obtain the spectrum of the example model in 
section [H 

References 

[1] D. Cremades, L. E. Ibanez, F. Marchesano, "Computing Yukawa couplings from mag- 
netized extra dimensions," JHEP 0405 (2004) 079. [hep-th/0404229]. 

[2] H. Abe, T. Kobayashi, H. Ohki, "Magnetized orbifold models," JHEP 0809 (2008) 043. 
[arXiv:0806.4748 [hep-th]]. 

[3] J. P. Conlon, A. Maharana, F. Quevedo, "Wave Functions and Yukawa Couplings in 
Local String Compactifications," JHEP 0809 (2008) 104. [arXiv:0807.0789 [hep-th]]. 

[4] J. J. Heckman, C. Vafa, "Flavor Hierarchy From F-theory," Nucl. Phys. B837 (2010) 
137-151. [arXiv:0811.2417 [hep-th]]. 

[5] H. Hayashi, T. Kawano, R. Tatar, T. Watari, "Codimension-3 Singularities and Yukawa 
Couplings in F-theory," Nucl. Phys. B823 (2009) 47-115. [arXiv:0901.4941 [hep-th]]. 



41 



[6] K. -S. Choi, T. Kobayashi, R. Maruyama, M. Murata, Y. Nakai, H. Ohki and M. Sakai, 
"E(6,7,8) Magnetized Extra Dimensional Models," Eur. Phys. J. C 67 (2010) 273 
[arXiv:0908.0395 [hep-pli]]. 

[7] A. Font, L. E. Ibanez, "Matter wave functions and Yukawa couplings in F-theory Grand 
Unification," JHEP 0909 (2009) 036. [arXiv:0907.4895 [hep-tli]]. 

[8] S. Cecotti, M. C. N. Cheng, J. J. Hecknian, C. Vafa, "Yukawa Couplings in F-theory 
and Non-Commutative Geometry," [arXiv:0910.0477 [hep-th]]. 

[9] J. P. Conlon, E. Palti, "Aspects of Flavour and Supersymmetry in F-theory GUTs," 
JHEP 1001 (2010) 029. [arXiv:0910.2413 [hep-th]]. 

[10] H. Hayashi, T. Kawano, Y. Tsuchiya, T. Watari, "Flavor Structure in F-theory Com- 
pactifications," JHEP 1008 (2010) 036. [arXiv:0910.2762 [hep-th]]. 

[11] L. Aparicio, A. Font, L. E. Ibanez, F. Marchesano, "Flux and Instanton Effects in Local 
F-theory Models and Hierarchical Fermion Masses," [arXiv:1104.2609 [hep-th]]. 

[12] V. K. Oikonomou, "F-theory Yukawa Couplings and Supersymmetric Quantum Mechan- 
ics," Nucl. Phys. B 856 (2012) 1 [arXiv:1107.0497 [hep-th]]. 

V. K. Oikonomou, "F-theory and the Witten Index," Nucl. Phys. B 850 (2011) 273 
[arXiv:1103.1289 [hep-th]]. 

[13] P. G. Camara, E. Dudas and E. Palti, "Massive wavefunctions, proton decay and FCNCs 
in local F-theory GUTs," JHEP 1112 (2011) 112 [arXiv:1110.2206 [hep-th]]. 

[14] G. Aldazabal, L. E. Ibanez, F. Quevedo and A. M. Uranga, "D-branes at singularities: 
A Bottom up approach to the string embedding of the standard model," JHEP 0008 
(2000) 002 [hep-th/0005067]. 

[15] J. P. Conlon, A. Maharana and F. Quevedo, "Towards Realistic String Vacua," JHEP 
0905, 109 (2009) [arXiv:0810.5660 [hep-th]]. 

[16] S. Krippendorf, M. J. Dolan, A. Maharana and F. Quevedo, "D-branes at Toric Singu- 
larities: Model Building, Yukawa Couplings and Flavour Physics," JHEP 1006 (2010) 
092 [arXiv:1002.1790 [hep-th]]. 

[17] C. P. Burgess, S. Krippendorf, A. Maharana and F. Quevedo, "Radiative Fermion Masses 
in Local D-Brane Models," JHEP 1105 (2011) 103 [arXiv:1102.1973 [hep-th]]. 

[18] M. J. Dolan, S. Krippendorf and F. Quevedo, "Towards a Systematic Construction 
of Realistic D-brane Models on a del Pezzo Singularity," JHEP 1110 (2011) 024 
[arXiv:1106.6039 [hep-th]]. 

[19] C. Vafa, "Evidence for F theory," Nucl. Phys. B469 (1996) 403-418. [hep-th/9602022]. 

[20] D. R. Morrison, C. Vafa, "Compactifications of F theory on Calabi-Yau threefolds. 1," 
Nucl. Phys. B473 (1996) 74-92. [hep-th/9602114]; 

D. R. Morrison, C. Vafa, "Compactifications of F theory on Calabi-Yau threefolds. 2.," 
Nucl. Phys. B476 (1996) 437-469. [hep-th/9603161]; 



42 



M. Bershadsky, K. A. Intriligator, S. Kachru, D. R. Morrison, V. Sadov, C. Vafa, "Geo- 
metric singularities and enhanced gauge symmetries," Nucl. Phys. B481 (1996) 215-252. 
[hep-th/9605200]; 

S. H. Katz, C. Vafa, "Matter from geometry," Nucl. Phys. B497 (1997) 146-154. [hep- 
th/9606086]. 

[21] R. Donagi, M. Wijnholt, "Model Building with F-Theory," [arXiv: 0802. 2969 [hep-th]]; 

[22] C. Beasley, J. J. Heckman, C. Vafa, "GUTs and Exceptional Branes in F-theory - I," 
JHEP 0901 (2009) 058. [arXiv:0802.3391 [hep-th]]; 

[23] C. Beasley, J. J. Heckman, C. Vafa, "GUTs and Exceptional Branes in F-theory - II: 
Experimental Predictions," JHEP 0901 (2009) 059. [arXiv:0806.0102 [hep-th]]; 

[24] R. Donagi, M. Wijnholt, "Breaking GUT Groups in F-Theory," [arXiv:0808.2223 [hep- 
th]]. 

[25] F. Denef, "Les Houches Lectures on Constructing String Vacua," [arXiv:0803.1194 [hep- 

th]]; 

J. J. Heckman, "Particle Physics Implications of F-theory," Submitted to: 
Ann.Rev.Nucl.Part.Sci.. [arXiv:1001.0577 [hep-th]]; 

T. Weigand, "Lectures on F-theory compactifications and model building," Class. Quant. 
Grav. 27 (2010) 214004. [arXiv: 1009.3497 [hep-th]]. 

[26] S. B. Giddings, S. Kachru and J. Polchinski, "Hierarchies from fluxes in string compact- 
ifications," Phys. Rev. D 66 (2002) 106006 [hep-th/0105097]. 

S. Kachru, R. Kallosh, A. D. Linde and S. P. Trivedi, "De Sitter vacua in string theory," 
Phys. Rev. D 68 (2003) 046005 [hep-th/0301240]. 

V. Balasubramanian, P. Berglund, J. P. Conlon and F. Quevedo, "Systematics of mod- 
uli stabilisation in Calabi-Yau flux compactifications," JHEP 0503 (2005) 007 [hep- 
th/0502058]. 

M. Cicoli, A. Maharana, F. Quevedo and C. P. Burgess, "De Sitter String Vacua from 
Dilaton-dependent Non-perturbative Effects," arXiv:1203.1750 [hep-th]. 

[27] T. W. Grimm and T. Weigand, "On Abelian Gauge Symmetries and Proton Decay in 
Global F-theory GUTs," Phys. Rev. D 82 (2010) 086009 [arXiv: 1006.0226 [hep-th]]. 

[28] T. W. Grimm, M. Kerstan, E. Palti and T. Weigand, "Massive Abelian Gauge Symme- 
tries and Fluxes in F-theory," JHEP 1112 (2011) 004 [arXiv:1107.3842 [hep-th]]. 

[29] A. P. Braun, A. Collinucci and R. Valandro, "G-flux in F-theory and algebraic cycles," 
Nucl. Phys. B 856 (2012) 129 [arXiv: 1107.5337 [hep-th]]. 

[30] J. Marsano and S. Schafer-Nameki, "Yukawas, G-flux, and Spectral Covers from Resolved 
Calabi-Yau's," JHEP 1111 (2011) 098 [arXiv: 1108. 1794 [hep-th]]. 

[31] S. Krause, C. Mayrhofer and T. Weigand, "G4 flux, chiral matter and singularity res- 
olution in F-theory compactifications," Nucl. Phys. B 858 (2012) 1 [arXiv: 1109.3454 
[hep-th]]. 



43 



[32] T. W. Grimm and H. Hayashi, "F-theory fluxes, Chirality and Chern-Simons theories," 
arXiv: 111 1.1232 [hep-th]. 

[33] S. Krause, C. Mayrhofer and T. Weigand, "Gauge Fluxes in F-theory and Type JIB 
Orientifolds," arXiv: 1202.3138 [hep-th]. 

[34] J. J. Heckman, A. Tavanfar and C. Vafa, "The Point of E(8) in F-theory GUTs," JHEP 
1008, 040 (2010) [arXiv:0906.0581 [hep-th]]. 

[35] S. Cecotti, C. Cordova, J. J. Heckman and C. Vafa, "T-Branes and Monodromy," JHEP 
1107 (2011) 030 [arXiv:1010.5780 [hep-th]]. 

[36] F. Marchesano and L. Martucci, "Non-perturbative effects on seven-brane Yukawa cou- 
plings," Phys. Rev. Lett. 104 (2010) 231601 [arXiv:0910.5496 [hep-th]]. 

[37] F. Marchesano, P. McGuirk and G. Shiu, "Open String Wavefunctions in Warped Com- 
pactifications," JHEP 0904 (2009) 095 [arXiv:0812.2247 [hep-th]]. 

[38] P. G. Camara and F. Marchesano, "Open string wavefunctions in flux compactifications," 
JHEP 0910 (2009) 017 [arXiv:0906.3033 [hep-th]]. 

[39] F. Marchesano, P. McGuirk and G. Shiu, "Chiral matter wavefunctions in warped com- 
pactifications," JHEP 1105 (2011) 090 [arXiv:1012.2759 [hep-th]]. 

[40] J. P. Conlon, "Gauge Threshold Corrections for Local String Models," JHEP 0904 
(2009) 059 [arXiv:0901.4350 [hep-th]]. 

J. P. Conlon and E. Palti, "Gauge Threshold Corrections for Local Orientifolds," JHEP 
0909 (2009) 019 [arXiv:0906.1920 [hep-th]]. 

J. P. Conlon and E. Palti, "On Gauge Threshold Corrections for Local HB/F-theory 
GUTs," Phys. Rev. D 80 (2009) 106004 [arXiv:0907.1362 [hep-th]]. 

[41] J. Marsano, N. Saulina and S. Schafer-Nameki, "Compact F-theory GUTs with U(l) 
(PQ)," JHEP 1004 (2010) 095 [arXiv:0912.0272 [hep-th]]. 

[42] J. C. Callaghan, S. F. King, G. K. Leontaris and G. G. Ross, "Towards a Realistic 
F-theory GUT," arXiv: 1109. 1399 [hep-ph]. 

[43] E. Dudas and E. Palti, "On hypercharge flux and exotics in F-theory GUTs," JHEP 
1009 (2010) 013 [arXiv: 1007. 1297 [hep-ph]]. 

[44] R. Donagi and M. Wijnholt, "Higgs Bundles and UV Completion in F-Theory," 
arXiv:0904.1218 [hep-th]. 

[45] J. Marsano, N. Saulina and S. Schafer-Nameki, "Monodromies, Fluxes, and Compact 
Three-Generation F-theory GUTs," JHEP 0908 (2009) 046 [arXiv:0906.4672 [hep-th]]. 

[46] E. Dudas and E. Palti, "Froggatt-Nielsen models from E(8) in F-theory GUTs," JHEP 
1001, 127 (2010) [arXiv:0912.0853 [hep-th]]. 

[47] S. F. King, G. K. Leontaris and G. G. Ross, "Family symmetries in F-theory GUTs," 
Nucl. Phys. B 838, 119 (2010) [arXiv:1005.1025 [hep-ph]]. 



44 



[48] G. K. Leontaris and G. G. Ross, "Yukawa couplings and fermion mass structure in 
F-theory GUTs," JHEP 1102, 108 (2011) [arXiv: 1009.6000 [hep-th]]. 

[49] P. Binetruy and E. Dudas, "Fermion mass hierarchies in low-energy supergravity and 
superstring models," Nucl. Phys. B 442 (1995) 21 [hep-ph/9411413]. 

[50] R. Donagi and M. Wijnholt, "Gluing Branes II: Flavour Physics and String Duality," 
arXiv:1112.4854 [hep-th]. 

[51] G. Ross and M. Serna, "Unification and fermion mass structure," Phys. Lett. B 664 
(2008) 97 [arXiv:0704.1248 [hep-ph]]. 

[52] J. Marsano, "Hypercharge Flux, Exotics, and Anomaly Cancellation in F-theory GUTs," 
Phys. Rev. Lett. 106 (2011) 081601. [arXiv:1011.2212 [hep-th]]. 

[53] M. J. Dolan, J. Marsano, N. Saulina, S. Schafer-Nameki, "F-theory GUTs with U(l) 
Symmetries: Generalities and Survey," [arXiv:1102.0290 [hep-th]]. 

[54] M. J. Dolan, J. Marsano, S. Schafer-Nameki, "Unification and Phenomenology of F- 
Theory GUTs with U(1)PQ," [arXiv: 1109.4958 [hep-ph]]. 

[55] L. E. Ibanez, C. Munoz and S. Rigolin, "Aspect of type I string phenomenology," Nucl. 
Phys. B 553 (1999) 43 [hep-ph/98 12397]. 

D. Lust, S. Reffert and S. Stieberger, "Flux-induced soft supersymmetry breaking in 
chiral type IIB orientifolds with D3 / D7-branes," Nucl. Phys. B 706 (2005) 3 [hep- 
th/0406092]. 

J. P. Conlon, D. Cremades and F. Quevedo, "Kahler potentials of chiral matter fields 
for Calabi-Yau string compactifications," JHEP 0701 (2007) 022 [hep-th/0609180]. 

[56] T. Kawano, Y. Tsuchiya and T. Watari, "A Note on Kahler Potential of Charged Matter 
in F-theory," Phys. Lett. B 709 (2012) 254 [arXiv:1112.2987 [hep-th]]. 

[57] L. B. Anderson, J. Gray, A. Lukas and E. Palti, "Two Hundred Heterotic Standard Mod- 
els on Smooth Calabi-Yau Threefolds," Phys. Rev. D 84 (2011) 106005 [arXiv: 1106.4804 
[hep-th]]. 

[58] L. B. Anderson, J. Gray, A. Lukas and E. Palti, "Heterotic Line Bundle Standard Mod- 
els," arXiv: 1202. 1757 [hep-th]. 

[59] E. Witten, "Supersymmetry and Morse theory," J. Diff. Geom. 17 (1982) 661. 



45 



